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Motivated by recent work on Hall viscosity, we derive from first principles the Kubo formulas for 
the stress-stress response function at zero wavevector that can be used to define the full complex 
frequency-dependent viscosity tensor, both with and without a uniform magnetic field. The formulas 
in the existing literature are frequently incomplete, incorrect, or lack a derivation; in particular, 
Hall viscosity is overlooked. Our approach begins from the response to a uniform external strain 
field, which is an active time-dependent coordinate transformation in d space dimensions. These 
transformations form the group GL(d, R) of invertible matrices, and the infinitesimal generators are 
called strain generators. These enable us to express the Kubo formula in different ways, related by 
Ward identities; some of these make contact with the adiabatic transport approach. The importance 
of retaining contact terms, analogous to the diamagnetic term in the familiar Kubo formula for 
conductivity, is emphasized. For Galilean-invariant systems, we derive a relation between the stress 
response tensor and the conductivity tensor that is valid at all frequencies and in both the presence 
and absence of a magnetic field. In the presence of a magnetic field and at low frequency, this yields 
a relation between the Hall viscosity, the q 2 part of the Hall conductivity, the inverse compressibility 
(suitably defined), and the diverging part of the shear viscosity (if any); this relation generalizes a 
result found recently. We show that the correct value of the Hall viscosity at zero frequency can be 
obtained (at least in the absence of low-frequency bulk and shear viscosity) by assuming that there 
is an orbital spin per particle that couples to a perturbing electromagnetic field as a magnetization 
per particle. We study several examples as checks on our formulation. We also present formulas 
for the stress response that directly generalize the Berry (adiabatic) curvature expressions for zero- 
frequency Hall conductivity or viscosity to the full tensors at all frequencies. 



I. INTRODUCTION 

There has been great interest recently in the viscosity 
of quantum fluids, coming from various directions. These 
directions include a conjectured lower bound on the ratio 
of the shear viscosity to entropy density of a fluid, from 
the AdS/CFT correspondence^; the properties of an in- 
teracting gas of fermions with interactions described by 
s-wave scattering at or near the unitarity limitS; and the 
so-called Hall viscosity, an antisymmetric part of the vis- 
cosity tensor analogous to Hall conductivity, which has 
been calculated for several gapped topological phases^—. 

Viscosity, whether in a solid or fluid, is essentially the 
set of transport coefficients describing the relaxation of 
a deviation of the momentum density from its value in 
(possibly only local) equilibrium. Hence it is necessary 
that momentum be conserved in order even to consider 
viscosity. If the momentum density at x at time t is 
g(x, i), the continuity equation for momentum is 



dt 



dfj.r^i^t) = 0, 



(1.1) 



(Greek indices /U, v, . . . , = 1, . . . , d refer to space com- 
ponents, d is the dimension of space, and repeated Greek 
indices are summed) . The stress tensor operator (x, t) 
plays a central role in calculating viscosity. The viscos- 
ity tensor in a fluid can be defined as the expectation 
of the stress due to a time- varying "strain" . Unlike in a 
solid, in a fluid with no external fields present, an intrin- 
sic local static strain tensor is not defined in full in any 



completely natural way, but we can define its trace as 
trit = 5(L d )/L d , where L d is the volume, for a uniform 
dilation, or locally using the average particle number den- 
sity n as tru = —5n/n. The local expected stress does 
not respond to a change in the shape of a box confining 
the fluid (which can be considered as an attempt to im- 
pose a static strain Uap), except that there is a response 
of the pressure to a change in volume; the pressure is the 
expectation of the trace of the stress tensor, divided by 
the dimension of space. However, time-dependent strain 
has an analog, which is the matrix of gradients of the 
velocity field v, the average velocity of the fluid: 



du a 



dt 



dx a 



The change in the average stress tensor (r Miy ) from its 
equilibrium value can be formally expanded in time 
derivatives, as 

5 (t m „) = —\f_ivapUap — -q^afj ^ + . . . , (1.2) 

where X^ap is the tensor of elastic moduli and rj^ap is 
the viscosity tensor. Here we may view this as holding 
between local quantities at the same position in space 
from a long- wavelength point of view, or to zeroth order 
in spatial derivatives. Then in a fluid, 



(1.3) 



that is, only the trace of r responds, and only to tru, 
corresponding to a change in volume or local density of 
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the fluid. The coefficient k 1 is the inverse compress- 
ibility, k- 1 = -L d (dP/d{L d )) N , the derivative of the 
pressure with respect to volume L d of the fluid, taken 
at fixed particle number N. So far, we did not use ro- 
tational invariance (isotropy) of the fluid; if we do, then 
the stress tensor is symmetric, and only the more familiar 
symmetrized rate of strain 

l_(dv« dv£\ (L4) 
2 \dxp dx a J 

enters, so that rj^ap is also symmetric under the ex- 
change a -R> /3. Another part of the expected stress is 
the momentum flux rrinv^v v (where m is the particle 
mass), which is very important in fluid mechanics (e.g. 
in obtaining the Navier-Stokes equations), but will not 
be encountered in this paper as we consider only linear 
response to strains. 

Without assuming rotation invariance, we can further 
distinguish some parts of the viscosity tensor. r\ can be 
divided into a symmetric and antisymmetric part under 
exchange of the first and the last pair of indices^: 

ii S = n S 

Vlivafi = ~VafJiJ.v (1-5) 

For the zero-frequency parts that we consider at the mo- 
ment, only the symmetric part of r] S contributes to dis- 
sipation. For a rotationally-invariant <i-dimensional sys- 
tem it has only two independent components, 

(1.6) 

with £ the bulk viscosity and r] sh the shear viscosity. r\ A 
can only be nonzero when time-reversal symmetry is bro- 
ken, and in rotationally-invariant two dimensional system 
has only one independent component, 

= V H (<W,x/3 - S Hi3e a u) , (1-7) 

with r) H the so-called Hall viscosity^!. For gapped quan- 
tum systems at zero temperature, it has been shown us- 
ing the adiabatic approach to quantum transport that 

r] H = \hns, (1.8) 

where n is the expected particle number density in the 
ground state, and s is minus the average orbital spin per 
particle^; see also Refs.Hand0- h is Planck's constant, 
which we usually set to 1. We should point out that in 
Ref. there is an unfortunate sign mistake in the defi- 
nition of the stress tensor in eq. (2.18) of that reference, 
which propagated through the paper, though the adia- 
batic curvature results are correct. This means that all 
viscosities calculated there should have the opposite sign; 
above we have also reversed the sign in the definition of 



the scalar rj , so as to retain the memorable formula, eq. 

(EHD. 

The purpose of this paper is to develop a variety of 
Kubo formulas from which the viscosity tensor, which 
in general is frequency dependent and complex, can be 
calculated in quantum fluids at zero or non-zero temper- 
atures. A motivation for doing so is to define the Hall 
viscosity from a Kubo formula, and thus make contact 
with traditional approaches. In addition, we consider the 
relation of the viscosity tensor to the conductivity tensor 
in Galilean-invariant systems, in which the latter tensor 
is the momentum-momentum density response; the rela- 
tion comes from the continuity equation, eq. (Il.l[) . and 
(along with some other relations in this paper) it can be 
called a Ward identity. This has been studied previously, 
or else frequently is used as the definition of the viscosity 
for calculation purposes; our goal is to do it allowing for 
the possibility of a Hall viscosity. Moreover, we will do 
the same in the presence of a magnetic field. In this case, 
total (kinetic) momentum is not time-independent, but 
precesses at the cyclotron frequency. We will show that 
we can nonetheless derive parallel formulas in this case. 
In particular, we recover a formula of Hoyos and Son£ 
that relates the Hall viscosity in a quantum Hall state to 
the order q 2 part of the Hall conductivity at wavevector 
q and the inverse compressibility kT 1 above (suitably de- 
fined in the presence of the magnetic field), and generalize 
it further. These results bring us closer to finding exper- 
imental techniques with which the Hall viscosity can be 
measured. We also study several example systems, to 
validate our approach to Kubo formulas for the viscos- 
ity tensor. Finally, we show, in a partially heuristic way, 
that the Hall viscosity can be rederived from macroscopic 
electrodynamics by using the relation with conductivity, 
and assuming that at low frequencies the system behaves 
as if there is an orbital spin — s per particle that acts as 
a contribution to magnetization density. 

The basic strategy of our work is to define the viscosity 
from the response of the stress to an external field. This 
external field enters (in some gauge choices) as a spatial 
metric, or if we do not assume rotational invariance, as 
what we will call a strain; it is generally assumed to be 
constant in space. (This strain field is external, and can 
have arbitrary time dependence. This does not contra- 
dict our statement that a static strain is not fully de- 
fined in a fluid, as that was for an intrinsic strain, which 
would be a property defined given any state of the fluid 
in no external field. But where the intrinsic strain is de- 
fined, i.e. for its static trace and for its time derivative, 
the stress responses to either of the two strains should 
agree.) It is well known (especially in gravitation the- 
ory and high-energy physics) that the stress tensor is the 
change in the Lagrangian or Hamiltonian with respect 
to the metric. We build on this to consider the under- 
lying response function, that gives the viscosity, as the 
next order response of this stress to the time-derivative 
of the strain field. This approach has the advantage of 
making contact with previous work in which the Hall vis- 
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cosity of some systems was calculated 3- ^ using the adia- 
batic transport technique. It closely parallels the case of 
conductivity, in which one considers the response of the 
current, which is the change in Lagrangian or Hamilto- 
nian with the vector potential, to an electric field, viewed 
as the time-derivative of the perturbing vector potential, 
setting that vector potential to zero at the end. Our 
approach uses the standard Hamiltonian (canonical) op- 
erator formalism, plus linear response. We do not make 
any assumptions that hydrodynamic behavior or local 
equilibrium holds, or use constitutive relations, beyond 
motivating the names for some parts of our expressions, 
such as compressibility and viscosity. We focus on quan- 
tum systems, but a similar discussion could be given in 
the classical setting. 

Early work on response-function type formulas for 
transport coefficients was done by Green, Kubo, Mori, 
Kadanoff and Martin, and Luttinger, starting in the 
1950s£~— . Green and Mori initially used somewhat phe- 
nomenological methods. "Mechanical" formulations as 
the response to an external field, similar to ours, were 
used by Kuboi2 for the electrical conductivity case, and 
by Luttinger 13 for thermal transport and viscosity, and 
for the most part the final formulas agree with those of 
Green and Mori. Kubo and Luttinger began by deriv- 
ing the retarded response function to an externally ap- 
plied field, and then transformed to expressions that are 
not as familiar today, involving an integration over an 
imaginary time variable, as well as one over real time. 
By contrast, most references from the last forty years 
follow Kadanoff and Martin^ 2 , in using only retarded re- 
sponse functions to obtain "Kubo-type" formulas, now 
usually called simply Kubo formulas. (We too use only 
this formalism.) For the particular case of viscosity, many 
authors including Luttinger 13 base the definition on the 
continuity equation, eq. (jl.ll) . and his "mechanical" for- 
mulation is a study of the response of the momentum 
density to an electric field (which couples to the number 
current, equal to the momentum density times the mass 
in a Galilean-invariant system of particles that all have 
the same mass and unit charge). This leads him to a cor- 
rect formula as a stress-stress response function which, 
however, he writes in terms of an additional integral over 
an imaginary-time interval. In many other papers, in- 
cluding the interesting recent Ref. [lj, the starting for- 
mula is the momentum-momentum response, which has 
to be expanded to order q 2 as the wavevector q tends 
to zero, to obtain the viscosity tensor. (We note that 
this approach cannot be employed to define viscosity of 
a finite system, which was defined using a strain in the 
adiabatic approach.) In some of these papers the stress- 
stress form is not used at all. But in most papers that do 
give a stress-stress form for the response function that 
yields the viscosity, the formulas frequently are incom- 
plete, lack a derivation, contain pitfalls for the unwary, 
or are simply incorrect. In particular, we are not aware 
of any derivations from an external strain field in the 
many-body literature. Hence, and because so much time 



has elapsed and notations have changed since the 1960s, 
we feel justified in revisiting these formal matters here. 

In order to explain the issues that need to be corrected 
in the stress-stress form, we will use the more familiar 
case of conductivity as an analogy. For the (complex) 
conductivity tensor at zero wavevector, it was shown by 
KuboiS that the approach outlined above gives rise to 
the Kubo formula 



in 



dte^ +t / d d x ([ju(x,t),j I/ (O,0)]) o . (1.9) 



<7/u/(w) 



Here j v (x.,t) is the current density operator (equal to 
g v jm in a Galilean-invariant system), u + = lu + ie, the 
limit e — > + is implicit, and the expectation is taken 
in the unperturbed ground state or statistical ensem- 
ble, denoted (. ..)o- We set the charge of the particles 
to 1 throughout this paper. This expression gives the 
response of the expectation of the current density to a 
uniform external electric field. (When mulitplied by m 2 , 
it is also the q — > limit of the momentum-momentum 
response mentioned above; this fact will not play any 
role just now). The first term is the so-called diamag- 
netic current, which (in the approach we are using at 
the moment) arises because the current density in the 
presence of a perturbing vector potential A v is actually 
j„(x,t) — A 1/ (x,t)n(x,t)/m (because the Hamiltonian is 
quadratic in A v ), where n(x, t) is the number density op- 
erator. The second term is the current-current retarded 
response function (though we usually reserve such termi- 
nology for the complete expression). We note that text- 
book derivations of linear response usually assume that 
the perturbing field appears only linearly in the Hamil- 
tonian (except possibly when considering conductivity), 
and then no contact terms appear. 

It would not be wise to drop the diamagnetic current 
term from the conductivity. Recall the Sokhotski-Plemelj 
formula (for u> real), 



— — = PP- - wt£(cj), 
u> + ie uj 



(1.10) 



where PP denotes the principal part (both terms become 
meaningful once substituted into an integral). Even if 
one wants to find only the real part of the conductivity, 
the diamagnetic term contributes a (^-function at zero 
frequency. In some cases, such as for fermions without 
impurity scattering (disorder), both in Fermi liquids and 
paired superfluids, the time-integral term vanishes and 
the diamagnetic term is the full response. At zero fre- 
quency, one would say that the conductivity is infinite, 
which is correct. 

For the real part, some authors instead use the formula 



Re<7 Ml/ (u;) 



Re f™dte^ +t Jd d x ([ > (x,t), > (0,0)]) 

(1.11) 
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in which the ttS(lj) coming from use of the Sokhotski- 
Plemelj formulas has been dropped, both in the dia- 
magnetic term, and in the time-integral term (the latter 
would contain the imaginary part of the time integral, in 
place of the real part here). This is correct when these 
(5-functions cancel, which does happen for the case with 
impurity scattering. In the simple Drude approximation, 
the complex conductivity tensor is 



<V,(w) 



in 



m(oj + i/ T imp)' 



(1.12) 



where l/ri mp is the impurity scattering rate; there is no 
S(uj) piece in the real part. In that case the real part of 
the time integral, which is kept, gives the correct broad- 
ened Drude peak, in the simplest approximation. But 
in the limit as the impurity scattering rate goes to zero, 
the correct result at zero frequency should be diverging, 
while use of the above formula gives a result that in- 
creases as the rate gets smaller, but jumps to zero when 
there is strictly no impurity scattering. One would like 
to criticize this behavior for not being continuous, how- 
ever if the value at zero scattering was infinity it would 
not be continuous either. Looking at the full function of 
ui, not only at u) = 0, we would like to say that the cor- 
rect result approaches a (^-function continuously as the 
scattering rate goes to zero. This makes sense only if 
we interpret Recr /i ^(w) as a distribution for any scatter- 
ing rate; then it is correct to say that as the rate goes 
to zero, it approaches a 5-function continuously (in the 
space of distributions). But this continuity with the limit 
of zero scattering is lost if the above form Re / °° dt . . . /ui 
is used. 

The use of eq. (|1.11[) also leads to difficulties with the 
Kramers-Kronig (KK) relations, that do not occur with 
the correct form, eq. (jl.9|) . even if the diamagnetic term 
is dropped (because that complex term obeys the KK 
relations). The difficulty can be remedied only by rein- 
stating the ttS(uj) times the imaginary part of the to = 
time-integral. However, the diamagnetic term cannot be 
recovered in this way, unless one has for example an ar- 
gument that the real part of a contains no (5-function at 
zero. 

The issues in the literature on Kubo stress-stress for- 
mulas for viscosity are very similar to these. Let us now 
give one of our forms (slightly simplified) for the response 
function from which we obtain the viscosity: 



1 



d\ a f: 



\=0' 



+i I dt I d d xe^ +t ([v,(x,i),r Q(3 (0,0)]) |>(1.13) 



The time-integral term is what would be expected for 
the stress-stress response. The first term, which is fre- 
quency independent except for the 1/uj + factor, is what 
we call (following the high-energy literature) a contact 
term. Without giving the full details here, A a( g is the ex- 
ternal strain, and varying the Hamiltonian with respect 



to A M „ gives (minus) the stress tensor r M „ (x) . The stress 
still depends on A, like the current above. Thus the re- 
sponse of the expectation of the stress to the strain con- 
tains the contact term, which is one more derivative with 
respect to \ a p ; it is directly analogous to the diamagnetic 
conductivity. By contrast, stress-stress response formu- 
las in the literature usually omit this term, an exception 
being the early work of Luttinger 13 , whose expression is 
equivalent to this, but is written in a way that may now 
seem obscure (the relation can be found in KuboiS). We 
will now give reasons why the contact term has to be 
retained. 

We pointed out earlier that in a homogeneous fluid, 
the expectation of the stress is 



(V,(x,£)) = P6^, 



(1.14) 



and that a static strain affects this at first order, the 
coefficient being the inverse compressibility. It follows 
that there will be a part of the response function, 



(1.15) 



as uj — > 0, and any remaining l/w + term has vanishing 
trace on the indices \iv and on a(3 (this result is valid 
at non-zero as well as at zero temperature). The coef- 
ficient kT 1 is equal to the zero frequency limit of the 
response function —iuj + Xtit*aa/d 2 , and there is a contri- 
bution from the contact term as well as from the time- 
integral term (we discuss this more formally in Sec. IIII.4I 
below). Thus dropping the contact term spoils this rela- 
tion. It is a feature of our approach that this term can be 
easily identified as the inverse compressibility, because it 
is the response to a static external strain. 

We then define the complex viscosity tensor to be 



(1.16) 



(We note that Luttinger— recognized the need to sub- 
tract such a term to obtain the viscosity tensor, however 
he derives it using hydrodynamic arguments and local 
equilibrium.) It follows from the preceding remarks that 
the bulk viscosity cannot diverge as u — > (see the def- 
initions above, for rotationally- invariant systems). The 
shear viscosity has a contribution from the contact term, 
and so the real part has a 5(uj) contribution. This may 
sometimes be cancelled by another from the time-integral 
term. In particular, this occurs in a gapped system at 
zero temperature, which should have no dissipative vis- 
cosity at uj — 0. It also occurs in an interacting Fermi 
liquid at positive temperature T, which has a finite shear 
viscosity at uj = that tends to infinity as the tempera- 
ture goes to zero. Thus, similarly to the case of the con- 
ductivity of a Fermi gas with impurities, the ^-function 
terms must cancel. Either of these cancelations is spoiled 
if the contact term is dropped. 
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Some authors, who consider fluids at positive tempera- 
tures, use an expression for the real part of the viscosity, 

Re / °° dt J d d x e^ +t ([t„„(0, t), r a[i (x, 0)]) 

u ' ( ' 

where u) in the denominator is real (ie is dropped). In 
particular, Kadanoff and Martin give a related form, 
which can be obtained from this using the fluctuation- 
dissipation theorem and letting ui — > 0. In this form, 
(5-functions S(ui) are all dropped, which removes the k _1 
term, and gives the viscosity correctly only if the remain- 
ing <5(u;)'s do cancel in the real part of our full expression 
above. As with the conductivity of a Fermi gas with im- 
purities, use of this form in an interacting Fermi liquid as 
T — » gives discontinuous behavior of the shear viscosity, 
which ought to be continuous when viewed as a distribu- 
tion. It also gives zero for a non-interacting Fermi gas, 
instead of infinity. While one may say that in either of 
these two limits, hydrodynamics is not well defined, it is 
preferable to have continuous behavior of our theories, 
and these are reasons to retain the contact term. 

A few authors go a step further than this, asserting 
that the complex viscosity is given by 

/ °° dt J d d x e- + ' ( [r M „ (0 , t) , r aP (x, 0)] ) Q 
w ' (L18) 

with 1/cj, not the more correct This is incorrect, 

as the right hand side usually has a divergence as uj — > 
in the trace part, which is a contribution to, but not equal 
to, k^ 1 . Further, if the 8(ui)'s are to cancel in the real 
part of the correct expression, it must also have al/u di- 
vergence in the imaginary part of the shear viscosity, that 
is cancelled by that due to the contact term. [The recent 
paper by Taylor and RanderiaAi correctly states that the 
contact term (in one particular form) contributes to the 
imaginary part of the viscosity] These latter forms also 
lead to difficulties with the KK relations, as in the con- 
ductivity case. There is sufficient information in these 
forms to recover a version (given by the last expression, 
but with ui replaced by lo + in the denominator) that sat- 
isfies the KK relations, but the contact term cannot be 
fully recovered in this way. 

A further reason to retain the contact term is that then 
the Ward identity relation with the conductivity tensor 
has a simple form, 

a v p(q,u)) = -^—8 ufj + q ^ qa +2 x^ai3(^)+0(q i ), (1.19) 

in zero magnetic field in a Galilean-invariant system. We 
see that there is a term iK^ 1 q^qp/(m 2 cu +3 ) on the right 
hand side, which is sometimes wrongly omitted (e.g. in 
Ref.0). 

Now wc will describe some of the highlights of our work 
in the present paper. We begin in Sec. [XT] with some tech- 
nical background used to set up the stress response ex- 
pressions. As mentioned above, we work systematically 



with the response to an external "strain" field. This can 
be introduced in two ways in particular. Again, an anal- 
ogy with conductivity may be helpful here. The conduc- 
tivity is the response of the number current density to 
an external electric field, and the field is taken to be uni- 
form in space, and have frequency to. The field can be 
represented in two ways in particular, which are simply 
two gauge choices. One is as a scalar potential, which de- 
pends linearly on position, because the field is uniform. 
The other is as a vector potential which is constant in 
space, and whose time derivative gives the electric field 
(and so is non-zero even when u> = 0). Both formula- 
tions appear in the literature; the latter is very common- 
place, while the first has the drawback that the linearly- 
varying scalar potential is not compatible with periodic 
boundary conditions, and the choice of suitable bound- 
ary conditions for a finite-size calculation becomes more 
problematic, though this is frequently ignored. On the 
other hand, in the first formulation the external electric 
field appears only linearly in the Hamiltonian, so there is 
no contact term in the response function, whereas in the 
formulation using the vector potential, the external field 
appears quadratically, and so there is a contact term — 
the diamagnetic current term, as discussed above. Last, 
we note that the transformation from one formulation to 
the other is a unitary transformation, simply implement- 
ing the change of gauge in the quantum theory. 

Similar alternatives appear in the case of stress re- 
sponse. It is known that stress is the derivative of the 
Hamiltonian with strain, so we may begin with a Hamil- 
tonian with a general, spatially-uniform strain that de- 
pends on time (it appears nonlinearly in the Hamilto- 
nian) . A time-independent strain can be eliminated from 
the Hamiltonian by a coordinate transformation, which 
for spatially-constant strain is a linear transformation. 
The two coordinate systems used to describe our system, 
which are related by this transformation, are referred to 
as the x and X variables, respectively. These coordinate 
transformations take the place of U(l) gauge transforma- 
tions; the relation can be understood if we realize that the 
conserved quantity corresponding here to particle num- 
ber in the conductivity case is the total momentum of the 
particles, and so in the generators of infinitesimal trans- 
formations, the x-dependence must be multiplied by the 
momentum of the particle on which it acts: 

^ = -^{4^} (1-20) 

i 

for zero magnetic field, which leads to coordinate trans- 
formations (for the conductivity problem, the corre- 
sponding operators are 'Yl li x l a ). These generators will 
be called strain generators; they play a central role in 
our approach. In fact, the time derivative of a strain 
generator gives the corresponding component of the in- 
tegrated stress tensor. This follows either from expand- 
ing the (Fourier-transformed) continuity equation to first 
order in wavevector, or by seeing that the time deriva- 
tive is the commutator with the Hamiltonian, which thus 
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produces the change in the Hamiltonian with a uniform 
strain, which we know is the (integrated) stress. 

As the strain generators generate linear transforma- 
tions of the coordinates (and inverse ones of the mo- 
menta, as the transformations are canonical), they obey 
the relations of the Lie algebra gl(d,M) of the group of 
invertible matrices GL(d, M) which describes (invertible) 
linear transformations. We further generalize this to the 
case with a magnetic field (in d = 2 dimensions). Pre- 
vious work 4,5 considered mainly transformations in the 
subgroup SL(d, R), consisting of shears and rotations but 
not dilations, because of the apparent difficulty of includ- 
ing dilations in the presence of a magnetic field. Here we 
overcome this difficulty using a technical trick: when the 
system is dilated, for consistency the magnetic field must 
be rescaled as an inverse length-squared. To do this, we 
promote the magnetic field strength to be a quantum 
variable B, for which the canonically conjugate momen- 
tum V does not appear in the Hamiltonian. Then oper- 
ator expressions for the strain generators can be written 
down as before, and the integrated stress tensor is again 
the commutator of these with the Hamiltonian. Finally, 
the states with which we work are assumed to have small 
width in B, so that the results correspond to a fixed mag- 
netic field as in conventional approaches. 

With this formalism in hand, we are ready to tackle 
stress response in Sec. IIIII One place to start is in the X 
variables, in which the time-derivative of the strain (but 
not the strain itself) appears linearly in the Hamiltonian, 
multiplied by the strain generator, similar to the conduc- 
tivity case when using the scalar potential to represent 
the external field; there is then no contact term in the re- 
sponse. This "stress-strain" form can then be integrated 
by parts on time, to obtain the "stress-stress" form which 
is similar to what was already discussed above, and in- 
cludes a contact term. (It can also be obtained as the 
response to a strain by working directly in the x vari- 
ables, similar to the vector-potential formulation in the 
conductivity case.) A different integration by parts pro- 
duces instead the "strain-strain" form, which is a strain- 
strain response function, and is useful in making contact 
with the adiabatic approach. These three forms have di- 
rect parallels with Kubo's Theorem 2 (there is a typo 
in the second line of the stated theorem in Ref. [HI the 
dot should be over the second occurrence of (f>BA, n °t 
the first). There are several variations on these results, 
including formulas for the response in the local stress 
instead of in the integrated stress, and the use of peri- 
odic boundary conditions in some of these formulations, 
or of a confining potential that is handled similarly to the 
magnetic field. We argue that the inverse compressibility 
must be subtracted from the foregoing response to ob- 
tain the viscosity tensor as described already above; in a 
magnetic field, this becomes the inverse "internal" com- 
pressibility, the partial derivative with area taken with 
the flux through the system fixed, not the field strength. 
We argue that the bulk viscosity is never diverging at 
zero frequency (unlike the shear viscosity, which may), 



and discuss the form of standard sum rules and positiv- 
ity constraints on the spectral density. We emphasize 
that these results are for all components of the viscosity 
tensor, including Hall viscosity. 

Sec. lIVl is devoted to the relation of the stress response, 
and hence viscosity, to the conductivity itself. The rela- 
tion is between the stress response at zero wavevector, 
and the second derivative of the conductivity with re- 
spect to wavevector. This has been discussed previously 
in zero magnetic field, and is sometimes used to define 
the viscosity. Our derivation includes a uniform magnetic 
field (with zero field as a special case) . This relation may 
then be studied at low frequencies, where if the bulk and 
shear viscosities are zero, there is a relation between the 
Hall viscosity, the second derivative of the Hall conduc- 
tivity, and the internal compressibility, which agrees with 
recent work by Hoyos and Son£. We find that the rela- 
tion still holds if these dissipative viscosity coefficients 
are non-zero but finite, and can be generalized to allow 
for diverging shear viscosity at zero frequency also. 

In Sec. |Vl we study several simple examples as checks 
on our formulation. These include the free (non- 
interacting) Fermi gas in zero magnetic field, which pos- 
sesses a diverging shear viscosity. Other examples are a 
non-interacting Fermi gas in non-zero magnetic field, for 
which we can find the full frequency-dependent response 
function, and fractional quantum Hall states, for which 
we recover the Hall viscosity in agreement with previ- 
ous work. Finally, we consider paired states of fermions 
for pairing functions with non-zero angular momentum, 
for which we recover the Hall viscosity in agreement with 
Refs.lj and[5|, for a mean-field model in which the Hamil- 
tonian is quadratic in field operators. In each case, we 
can verify the relation with conductivity, using results of 
some previous authors (who did not consider viscosity). 

Finally, in Sec. I VII we show that using macroscopic 
electrodynamics, together with the assumption (in the 
spirit of an effective theory) that an external perturb- 
ing magnetic field couples to the orbital spin — s of each 
particle as a magnetization, we again recover the Hall 
viscosity result, eq. (|1.8p . 

In Appendix we describe a formalism for including 
a confining potential to make a finite-volume system in 
an infinite space, and show that it works in the example 
of free fermions in a harmonic potential. In Appendix |B1 
we show that the three forms of Kubo formula for stress 
response can be written down in periodic boundary con- 
ditions, using derivatives of the ground state instead of 
strain generators. This makes contact with the adiabatic 
approach to quantum transport, in which Hall conduc- 
tivity and Hall viscosity are written as Berry curvatures 
(and then as a Chern number in the case of Hall conduc- 
tivity) and generalizes it to non-zero frequency, and to 
gapless systems. In Appendix [Cl we show that some as- 
sumptions of time-translation invariance that were used 
at some points in the derivations are justified. 
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II. STRESS TENSOR AND STRAIN 
GENERATORS 

In this section, we provide background needed for the 
viscosity linear response calculation. We discuss Hamil- 
tonians with a strain that could be static, and others 
(related by a coordinate transformation, implemented by 
a unitary transformation) that contain only the time- 
derivative of the strain, times a strain generator. The 
basic relation giving the integral of the stress tensor as 
the time derivative of the strain generator is obtained. 
We consider separately the two cases of zero magnetic 
field (in Sec. III.l[ >. which is mostly straightforward, and 
non-zero magnetic field (in Sec. III.2[) . which is less so. 



II. 1. Zero magnetic field 

We begin by considering a Hamiltonian H a for a system 
of interacting particles in infinite d-dimensional space: 



-y 

2m 4^ 



PuPu 



x J ), (2.1) 



where x* are the coordinates of the ith particle (i = 1, 
. . . , N) , and p 1 are their momenta, with 



= iS,,,,d 



(2.2) 



(Here and in what follows, we use Roman letters i, j, . . . , 
for particle indices.) We could be more general by intro- 
ducing an anisotropic mass in the kinetic term; this mod- 
ification is simple to make and will not be done explicitly. 
In general we do not assume the interaction potential V 
is rotationally invariant. Then we introduce a spatially- 
uniform strain A, that is a linear transformation in the 
coordinates (with matrix A M „), and the corresponding 
inverse transformation of the conjugate momenta^. We 
then allow A to depend on time t. This transformation 
is viewed actively, as changing the Hamiltonian of the 
system, relative to the same variables x l , p l . The Hamil- 
tonian for this system is then 



H A (t) 



where 



2m 



E^ + |E y ( AT w( xl - xJ ))' 



«r(*)=A-i(*)A-i(t). 

We also define the matrix A p „ by 



A = e> 



(2.3) 



(2.4) 



(2.5) 



as matrices. Here g^ is the metric, and g^ v is the inverse 
metric. If the Hamiltonian is rotationally invariant, then 
the strain enters only through the metric g^ u and its 



inverse. We shall work primarily in units where h = 
e = 1, unless otherwise stated. 

The strain transformation is canonical (and in fact 
does not mix coordinates and momenta), so it can be 
implemented by a unitary transformation S(t) for each 
i: 



(2.6) 



We can parametrize the strain transformation S(t) in 
terms of generators J^ v such that 

S(t) = exp(-iA^(t) V). (2.7) 

Inserting this into Eq. (|2.6[) we find that the self-adjoint 
strain generators J must obey 

It follows that the strain generators must also satisfy 

i [JfiUi Jafs] = b^fiJav — SvaJufS, (2.9) 

which are the commutation relations of the Lie algebra 
gl(d, R), associated with the group GL(d, R) of linear co- 
ordinate transformations. In particular, the antisymmet- 
ric part of J^ v is simply —1/2 times the angular momen- 
tum operator; for example in three dimensions, 

L a = — toiiv J^v (2-10) 
We can satisfy Eqs. (|2.8p and (|2.9[) by choosing^ 

i 

With H A (t) as the starting point, we now use a similar 
strain transformation but with a different point of view. 
We make a time-dependent change of variables from x 1 , 
p l (which we term the x variables) to X 1 = A T (t)x I and 
P 1 = A _1 (t)p 4 (which we term the X variables), which 
is again a canonical transformation implemented by the 
same S is defined as above: A T (t)x. 1 = S(t)x l 5 _1 (i), 
and similarly for p\ In the Heisenberg picture, the time- 
dependence of operators is determined by the Heisenberg 
equation of motion, for example for an operator A in the 
x variables 



i— = [A,H A (t)]+i- 



df 



Of 



(2.12) 



where dA/dt means the derivative acting on the explicit 
time dependence of A. This equation of motion requires 
a choice of the canonical variables used (here x l , p z ), 
which are viewed as having no explicit time dependence. 
But then due to the time-dependence of the change of 
variables, when the X variables are viewed as having no 
explicit time dependence, the resulting equations of mo- 
tion will not agree with simple changing variable unless 
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an additional term is included in the Hamiltonian. For 
an operator B, using the X variables, one has 



dB . . .8B 
l — = [ B) H}+i— ) 



dt 



at 



where 



(X) 



H = H, 



#1 = -iS- 1 — 



Hi 

as 
at 



(2.13) 

(2.14) 
(2.15) 



x,p 



in which Hq X> = H\ (t) is viewed as a function of X 1 and 

P 1 , and so is the same functional form as above, 
but with x 1 and p l replaced by X 1 and P l . In the partial 
(explicit) time derivative of S, x l and p 1 are to be viewed 
as time-independent. However, one can show that 



dS 

di 



— S-\t). (2.16) 



x,p 



S can be expressed in terms of strain generators for the 
X variables, 



t(x) 



(2.17) 



which are related to by a linear transformation. The 
change of canonical variables is similar to that which 
occurs when passing to a rotating frame in mechanics, 
which in fact is a particular case of the above derivation. 

As we will be interested only in linear response, we can 
expand the perturbation 
find that in the X variables the term Hi is 



iS x ^ to first order in A to 



rr _ 9 V r 



(2.18) 



(This form of a perturbation to a rate of strain is fairly 
common in the literature, see e.g. Ref. [T3-) At the same 
time, the distinction of x and X can be dropped, and 
we usually use x variables to lighten notation. We have 
thus mapped the Hamiltonian for a system with a time 
varying metric to that of a system in a fixed Euclidean 
metric, with a perturbation generating time-dependent 
deformations. The viscosity can now be computed as the 
response of an appropriately defined stress tensor to this 
strain perturbation. This is analogous to the computa- 
tion of the conductivity , where one can consider ei- 
ther the response of the current to a time- varying vector 
potential, or — after an appropriate gauge transforma- 
tion — the response of the current to a scalar potential. 

To define the stress tensor, let us start with the con- 
tinuity equation for momentum density g(x), defined in 
zero magnetic field and in the absence of strain (and so 



using x variables and Hamiltonian Hq ) by 



;(x) = i£{pV(x-x*)} 



(2.19) 



In the absence of other external forces, the continuity 
equation reads 



^^+«)(x,0 = 0, 



(2.20) 



where t$ (x) is the stress tensor operator in the absence 
of strain, and <9 M = d/dx^. Expressions for t$) for the 

Hamiltonian can be found in Refs. [TH and [TH We 
note that the continuity equation does not completely 
determine rju) , because any operator with vanishing di- 
vergence (such as a curl of something) could be added to 
t without violating the equation. This should not affect 
physical results, and in particular will not be an issue 
when the stress tensor is integrated over all space, as it 
often will be when the viscosity is calculated. 

Because we are interested in the long-wavelength be- 
havior of the stress tensor, we will examine Eq. (|2.20[) in 
Fourier space. To leading order in wavevector q, we have 



ig„7#>(q = 0)=0. 



We define T$) = T$(q = 0) or equivalently T$) = 
J d d x Tfw (x), the integrated stress tensor. Then, because 




Hq is translationally invariant, 



T^) = -^=-i[H ,J llv \. 



dt 



(2.21) 



This relation between the stress and the strain generators 
is a central result of this section. (As an example, for non- 
interacting particles, Tfi/ = Y^iP i uPl/' m > the momentum 
flux.) Since the antisymmetric part of J^ v is the angular 
momentum, it is clear that the integrated stress tensor 
Tf£) is symmetric (that is, the antisymmetric part is zero) 
when the Hamiltonian is rotationally invariant. We can 
also view the result in another way: referring to eq. (12.61) , 
we see that 



T, 



(o) _ 



dH A 



ox 



fJLV 



(2.22) 



A=I 



This is consistent with the idea that the stress tensor can 
be obtained by varying the Hamiltonian with respect to 
the metric, but here is slightly more general as we do 
not assume rotational invariance. In the rotationally- 
invariant case, we can also write 



T (o) 



, dH A 
J dg^ ' 



(2.23) 



which again is clearly symmetric, because is so (we 
note that the corresponding expression in Ref. |5j has the 
wrong sign). 

Strictly speaking, because the Hamiltonian H con- 
tains no confining potential, for repulsive or for no inter- 
actions, there is no normalizable ground state. We will 
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nonetheless go ahead somewhat informally (as in many 
papers in the literature), and evaluate the response in a 
normalizable state in which the particles occupy a region 
(or "box") of volume L d , and we will treat the state as 
if it were an eigenstate. Such expectation values will be 
written as (. . .)o. (A similar approach can be used for 
non-zero temperature, and most of the following is also 
valid for that case.) In a large system (i.e. as L and N 
become large with the density remaining fixed — we re- 
fer to this as the thermodynamic limit) this procedure 
may possibly be justified over intermediate time scales. 
In any case, our results will be seen to make sense. More 
formally, and completely generally, we can add to H\(t) 
an explicit "confining" potential term U, 



[/ = 5>(x'). 



(2.24) 



The single-particle potential u(x) has no A dependence, 
and can be used to represent a box (say, by using a hard- 
wall potential) the shape of which is fixed in x variables 
even as A varies. Then in the X variables, the poten- 
tial becomes U = ^ u(A T X), so the shape of the box 
varies with A (these conventions agree with Ref. d, ex- 
cept that there a periodic boundary condition was used 
instead, and the shape was chosen to be a square in the 
x variables). The potential modifies the continuity equa- 
tion for momentum density by the inclusion of an ex- 
ternal force due to the potential, while the stress tensor 
remains unchanged. We can also extend our formalism 
to include the potential, in such a way that an analog of 
cq. (|2.21[) holds, with H + U in place of H , by modi- 
fying the strain generators. In this way, we can justify 
all the later results that we present here less formally by 
ignoring the potential that confines the system to a finite 
volume. This is sketched in Appendix 1X1 

Returning to the original approach without a confining 
potential, at zeroth order in strain, we can consider the 
expectation of the stress. Then use of eq. (|2.22p and the 
Fcynman-Hellman theorem gives 



/o 



dE{\) 



A=// N 



(2.25) 



where E(X) = (Ha)q is the energy eigenvalue of the A- 
dependent ground state of H\ (for time- independent A) , 
and we choose to emphasize that particle number N is 
held fixed. By evaluating E in the X variables, in which 
a time-independent A enters only through the shape and 
size of the box, we realize that when the state is a ho- 
mogeneous fluid, the ground state energy depends on the 
volume of the box but not on its shape (in leading order 
in the thermodynamic limit). As —dE/d(L d ) — P, the 
pressure, we have therefore 



(T$) =PL% U , 



(2.26) 



(to leading order) as would be expected from the stan- 
dard result, (r^(x))o = PS^. For non-zero tempera- 
ture, we obtain the same result, as the derivative of the 



expectation of energy is taken with the probabilities held 
fixed, in order to use the Feynman-Hellman theorem un- 
der the thermal average. This implies that the derivative 
is taken with the entropy fixed, in agreement with the 
thermodynamic definition of pressure. (A similar method 
is used to find an expression for the pressure in Ref. [14J.) 
There are certain subtleties about the argument for this 
result (for either zero or non-zero temperature), which we 
will briefly discuss later in Sec. IIII.4l however the result is 
still valid. We emphasize that we did not use rotational 
invariance to obtain it. 

We also point out that if we use eq. (|2.2ip and then 
argue that the ground state is an eigenstate of H , we 
will conclude that all components of (T$))o vanish. This 
is incorrect in general because in the absence of the po- 
tential [/, the normalizable state that we use is not an 
eigenstate of Ho for repulsive or no interactions, while 
if the potential U is included, there is a normalizable 
ground (energy eigen-) state of H + U, but again the 
argument is blocked. 

In order to derive linear response of the stress to a time- 
dependent strain, we will need the stress to next order 
in the strain. For this it is convenient to notice that 
H\ is translationally invariant, and so the sum of p l is 
conserved. This means that the density of p momentum 
obeys a continuity equation, even in the presence of the 
time-dependent strain A(i): 



dt 



(2.27) 



Here t^ v depends on A; to illustrate the form of this, note 
that the momentum flux term part of integrated over 
space is 



(2.28) 



where the omitted terms refer to interactions. Following 
the same derivation as before, we see that 



(2.29) 



Then expanding to order A, we have 



T, 



r(x) 



(2.30) 

However, when we calculate response in the X variables, 
we will use the stress tensor in those variables, which is 
defined by transformation of tensors, so 



In terms of the X variables, this has the form T, 



(2.31) 

(x) _ 

Si PliPll m + • • ■ (where the omitted terms are from in- 
teractions), and coincides with T$ with X 1 and P l in 
place of x l and p\ In the x variables, it can be expanded 
as 

^ - -*[^ X) , VI + Kf3[[H { *\ V], J a p] + 0(A 2 ), 

(2.32) 
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in which the order in the double commutator is reversed, 
compared with T^ v . It is this integrated stress T^ v that 
we believe constitutes a natural starting point for the 
calculation of viscosity. However, we will see that in a 

homogeneous fluid, the difference in results from using 
fx) 

either T M „ or T^ v is negligible in the thermodynamic 
limit. 



II. 2. Nonzero magnetic field in two dimensions 

We now turn to the problem of generalizing the preced- 
ing set-up in the presence of an external magnetic field 
B] we concentrate on two dimensions. Most of the work 
is in finding the strain generators for this case; we will do 
so by considering separately the cases of pure shear defor- 
mations (det A(t) = 1) and pure dilations (A M „(i) oc 8^ v ). 
The unstrained Hamiltonian can be taken to be 



x?), (2.33) 



where 



7r l =p 4 -A(x*), 

[7rj,7r*] =iB5 3k e^. (2.34) 

We usually assume that the interaction V is independent 
of the magnetic field. Following the procedure of Sec- 
tion |ILT] above, we seek strain generators J M „ satisfying 



^ \jlJ>v ) x "a\ — 3av%'ji 



(2.35) 



First, let us consider the case of pure shear deforma- 
tions. The Hamiltonian for the system in the presence of 
the time- varying strain with det A = 1 is given by 

H *® = h E ^(tKK + lE v ( aT - xJ )) ■ 



(2.36) 

The condition det A = 1 implies that trA = 0. If we 
attempt naively to generalize the generators from zero 
magnetic field by taking 



tfw — E 9 { 2 7" 7r ^} 



(2.37) 



we find that 



/sh j 
u [list "a 



If we define 



S = cxp (-itr(A T J sh )) 



then these relations imply that the coordinates x 1 trans- 
form correctly, while the momenta tt 1 transform as 



0(X 2 



But note that if A Mi , is traceless, then the curl of the extra 
term containing B is 0, and so it is just a A dependent 
gauge transformation. Thus, we see that J generates 
the desired strain transformation along with the gauge 
transformation 

A^ y A^ BXaj^Ci/^Xa. 

An integration of this gauge term shows that, since A is 
traceless, the gauge transformation is given by 

In order to remove this unwanted gauge transformation, 
as well as to make this J^ v traceless, we define 



Jf v = jfu- Jtr(J sh )<V + f E e 



(2.38) 



E 



B 



A short calculation shows that defined in this way re- 
produces the traceless part of the transformations (|2.35|) , 
and therefore it is the desired traceless strain operator. 

Next, we consider the case of a pure dilation. We will 
soon see that we must define strains of the system so that 
they rescale the magnetic field in such a way that the 
magnetic flux $ = L 2 B through the system stays fixed 
while its shape is strained. With a fixed particle number 
N, and defining the filling factor for the region of area L 2 
occupied by the particles as v = 2ttN/ (BL 2 ) (as usual), 
this has the effect that we consider deformations at fixed 
filling factor v. For a dilation, we have 



A,, 



1 



tr(A)5 M „. 



In accordance with Eq. (|2.35p . we seek a dilation gener- 
ator K satisfying 

i[K,^] = t£. (2.39) 
Writing ir' 1 ^ = p^ — (x* ) , we see that these imply that 

i [K,pl] = P i li) 
i[K,A^)]=A^). (2.40) 

At first glance, it appears that we are in a quandary — 
the first and third of these equations are inconsistent, un- 
less the magnetic field strength B also transforms under 
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the action of K. To accomplish this, we promote the field 
strength B to a dynamical variable, represented by an op- 
erator B, whose eigenvalues are values B. This enlarges 
the Hilbert space of the system to include states with dif- 
ferent magnetic fields; we will continue to consider only 
B > 0. At this point, we are naturally motivated to 
introduce the momentum V conjugate to 23, such that 



[B,P] 



(2.41) 



B and V commute with x l and p* for all i. Note that the 
vector potential is now a function of B 7 but is independent 
of V . Wc then have 



A^\B), 



(2.42) 



and 



[7rj,7r*] = iBS jk e^. (2.43) 
For consistency with eq. (|2.39[) . we must also have 

i[K,B]=2B. (2.44) 
The Hamiltonian (|2.36j) now becomes 

+ i^(A T (t)(x'-x% 



(2.45) 



where we have expressed 7r in terms of the canonical mo- 
mentum and vector potential in order to make explicit 
the dependence of the Hamiltonian on B. (Again, an 
anisotropic mass can be introduced in the kinetic terms 
if desired.) Since H is independent of V, states with 
given eigenvalues B of B retain those values for all times. 
Eigenstates of B are not normalizable, however, instead 
we can use normalized packets with very small width in 
B to calculate expectations in linear response. 

A gauge-invariant dilation generator satisfying the 
commutation relations (|2.39[) and (|2.44[) is then given by 



K = — 



+ {B,S({x < },B)} 



(2.46) 



where we have introduced the "kinetic momentum" 
S({x*},Z?) = V - J2 t A x \ B ) conjugate to B, and 
{x 1 } = {x l : i = 1, . . . , N} is the set of all x"s. Under a 
gauge transformation generated by some scalar function 
(Mx, B) we have 



A M (x, B) -> A„ + M 0(x, B), 
A(x, B) -> A + <9 B 0(x, B). 



(2.47) 



To complete our definition of K, we specify that if 
the gauge choice is the symmetric gauge for all B, 
that is A M (x) = —\Be^ v x v (which is preserved by x- 
independent gauge transformations), then A is a function 
of B only, independent of x. Then 



[<,{/3,S({x fc },/3)}] 



Br 3 



(2.48) 



and this result of course is gauge covariant. 

Putting it all together, we have thus shown that 



7 _ /sh 



=J2 H-{ x i><} +Be ^ x i<) 



-V{S,S({x 4 },6)} 



(2.49) 



gives the strain generators for two dimensional systems in 
a magnetic field, which satisfy eqs. (|2.35|) . (|2.44[) . These 
generators also satisfy the g[(2,R) commutation rela- 
tions (|2.9p . and the antisymmetric part is —1/2 times 
the gauge-invariant rotation generator ( "angular momen- 
tum"). We note also that we consider only the range 
B > 0. It can then be shown that the operator {B, S} 
is a bona-fide self-adjoint operator, which is not the case 
for the operator S alone, on this range. 

As in Section III. 11 we can now apply the canonical 
transformation S = cxp(— iX^ J^ v ) to transform to the 
Hamiltonian in X variables, in which the corresponding 
kinetic momenta are II 1 (the variables B and V in X vari- 
ables should be distinguished from those in the x vari- 
ables also, but we will not introduce additional notation 
for this; at this point, it should be clear from the con- 
text). Wc find, as in the zero magnetic field case, that 
the system with time- varying strain is equivalent to one 
with Hamiltonian 



dt 



(2.50) 



to first order in A. 

As above, we can obtain an expression for the stress 
tensor in the presence of a magnetic field by considering 
the continuity equation for the kinetic momentum den- 
sity (with time-dependence obtained from , which is 
now generalized to include B), which reads 



^+^)(x)^e raJa (x), 
at H m 

with the kinetic momentum density given by 



g( x ) 



(2.51) 



(2.52) 



As in the previous section, we Fourier transform this 
equation, and to first order in wavevector q we find 
(again, T$ =Jd 2 xT$( x )) 







«7m 
2 



i(E(<- ! f i H'<})]+^) 



m 



(2.53) 
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Now, the first term on either side cancels since only the 
Lorentz force breaks conservation of total kinetic momen- 
tum. Thus 

= ^E (| K><} - • (2-54) 



We see that the term in the time derivative matches the 
first term in Eq. (|2.49l) for the strain generator. Then we 
have 



d_ 

di 



(v + ^EKx}) = 



B 



7T- E ( e » a (^a + X l n a) + V^a^a) 



2m 
B 



2^ Q EKX} 



Thus finally 



T (o) 

[11/ 



dt 



rr( x ) t 
n ' J uv 



(2.55) 



(2.56) 



just as in Eq. (|2.21|) above. The definitions for T M „ and 

for Tj^n and their expansions to order A in x variables 
are the same as in zero magnetic field. 

Again, at zeroth order in the strain, we can express 
the expectation of the stress in terms of thermodynamic 
properties, if we assume the state is a normalizable 
ground (and eigen-) state of Ho for each value of B, and 
as a function of B is concentrated near a value B. In 
the case of a two-dimensional system with a magnetic 
field, the field itself provides confinement, so normaliz- 
able eigenstates exist for given B, and we may consider 
a disk of fluid that covers an area L 2 . Then we expect 
that the average total stress (T$))o can be decomposed 
into a pressure term and a magnetization termi^ as 



MB), 



(2.57) 



where P = — (dE/d(L 2 ))N : B is the pressure, and M — 
— (dE/dB)]six 2 is the total magnetization. We define the 
internal pressure^ as 



Pint = P 



MB 



by defining the internal pressure, or equivalently by tak- 
ing the derivative with the flux through the system held 
fixed. 

In our framework, derivatives with respect to strain are 
taken in exactly that way. If we calculate the expectation 
of the stress using eq. (|2.56[) . then we must be careful 
to recall that the state is a wavepacket in B, and so 
not an energy eigenstate of Hq, because the energy for 
given B generally depends on B. Now when calculating 
the expectation of an operator that may depend on B, 
but not on V (i.e. does not contain V), we can take the 
expectation in the Hq eigenstate for each B, and then 
average the result over B using the dependence of the 
wavefunction on B. For the traceless (or shear) part of 
Juv, V does not appear, and so the fact that the state is 
an eigenstate of Ho for each B can be used to conclude 
that the expectation of T$ is zero. This cannot be done 
for the trace, so the result can be non-zero. On the other 
hand, the stress, including its trace, is itself independent 
of V , and so the result is the average over a small range of 
B of the result for each B, and the latter can be related 
to the derivative of energy with the strain (holding v 
and N fixed) using the Feynman-Hellman theorem just 
as in the zero magnetic field case. (This is done in x 
variables, using H\, and we emphasize that then x l and B 
are viewed as fixed when taking derivatives with respect 
to strain A.) Then by averaging over a sufficiently small 
range of B, the result for the trace is just the internal 
pressure. Then the full result is 



hit 



(L)L 2 



(2.60) 



(2.58) 



and this result can be considered as exact, rather than 
just as the order L 2 part as the size goes to infinity, if 
Pint(L) is defined in this way at finite L, but tends to the 
thermodynamic P m t as L — > oo. For non-zero tempera- 
ture, we should include all energy eigenstates, weighted 
by their Gibbs weight, but then a similar difficulty as 
in the zero magnetic field case in infinite size reappears: 
most states are spread over arbitrarily large volumes. We 
may deal with this in the same way as in Sec. 111.11 or by 
including a confining potential as in Appendix [X] The 
result takes the same form as in eq. (|2.60p . 



III. KUBO FORMULAS FOR VISCOSITY 



or as 



Pi, 



/ dE 



(2.59) 



where the partial derivative is at fixed TV and fixed filling 
factor, as we specified before. As pointed out by Cooper 
et. alJ£, in a homogeneous fluid, the usual pressure P is 
the change in energy under a change in the size of the 
box (at fixed B), and so includes a contribution from 
the Lorentz force acting on the boundary current that 
is related to the magnetization. This part is removed 



In this section, to calculate the viscosity, we first com- 
pute the retarded response function X^ a p of the inte- 
grated stress tensor T^ v to the perturbation H\ in Sec. 
IIII. II Initially, we consider separately the zero magnetic 
field case in any dimension and the case of nonzero mag- 
netic field in two dimensions. In Sec. IIII. 21 we give the 
response in the form of an intensive response function 
X, which can be defined in periodic boundary conditions 
also. Then in Sec. IIII.3l we relate X and x to the viscosity 
tensor. We show in Sec. IIII.41 that a leading part of x a t 
low frequency is the inverse compressibility, and finally 
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in Sec. IIII.5I discuss sum rules and positivity properties 
of X and x- 



III.l. The response function from strain generators 

We work in the X variables, and calculate the response 
fx) 

of the integrated stress T^ u ' to the perturbation Hi, us- 
ing linear response theory. Dropping the superscript (X), 
the change in the stress tensor to first order is given by 

(T^>(t)-<T^) = - f dt'X^{t-t')^P-, 

(3.1) 

where the retarded response function X is given by 
X^ aP {t) = - lim i0(t) ([Tp„(*), J Q ^(0)]) o e- et . (3.2) 

These expressions are to be evaluated with A — 0, so 
we may now drop the distinction between X and x, and 
corresponding superscripts. The time evolution here is 
taken with respect to Hq, and the expectation is taken 
in the unperturbed ground state of Hq (in zero magnetic 
field, it is again subject to the same caveats as in sec- 
tion III. 1 1 which are addressed further in Appendix [A} . 
The exponential damping ensures that the system was 
unperturbed infinitely far in the past. Fourier transform- 
ing Eq. (|3.2[) . we find that in the frequency domain 



X 



lim i 

>°+ jo 



dte^ +t ([T^(t),J a/3 (0)]) ( 



(3.3) 

where oj + = oj + ie (the e — > + limit will be left implicit 
from here on). We call this the stress-strain form of the 
response function. It corresponds directly to the response 
of the stress to an applied spatially-uniform rate of strain 
d\ a /3 jdt. No rotational invariance has been assumed in 
the derivation. 

Using the relation (|2.21l) . we can express Eq. (|3.3[) in 
two additional equivalent forms. The second form of the 
Kubo formula is the stress-stress form 



1 



PV(0),J a/3 (0)]> 



+ jf° dte^ +t ([T^(t),T aP (0)])o) ■ (3-4) 

(In obtaining this, we had to use time-translation invari- 
ance of the correlation function in eq. (|3.3p to shift the 
time dependence onto the operator J a p, and then back 
after using the identity (|2.2ip and integrating by parts 
on t.) The time-integral term is what one might have 
expected for the response function, as the deformation 
of shape couples directly to the stress tensor through 
the metric. In the additional contact term (the equal- 
time commutator term), the coefficient of 1/uj + is purely 
imaginary, because it is the expectation of a commutator 



of self-adjoint operators. The complete expression is di- 
rectly analogous to the standard Kubo formula for con- 
ductivity in terms of the current-current response, and 
the contact term in eq. Q3.4p is analogous to the diamag- 
netic conductivity. The latter will be discussed further 
in Sec. HTTT21 below. 

Lastly, we have the strain-strain form of the response 
function 

Xfu, a ff(u) = - i ([J M v(0), J a p(0)]) 

+ oj+ / dte w+t ([J^(t), Jafi(0)]) . (3.5) 
Jo 

(In this case, the identity (|2.21l) and integration by parts 
was used on the operator T, IU in eq. 



and time- 
translation invariance was not invoked.) This form of 
the response function is closely connected with the adi- 
abatic formalism for viscosity^—. For systems with non- 
degenerate ground states and an energy gap, the contact 
term (the equal-time strain-strain commutator term) is 
the adiabatic curvature associated with deformation of 
the metric, and gives the full response asw^O. Note 
this part is real and is manifestly antisymmetric under 
exchanging the pair \iv with aft. It gives the simplest 
way to see that the Hall viscosity is connected with the 
orbital spin density in these cases^. See also Appendix 

m 

A similar analysis applies in the magnetic field case. 
Because the equations of motion Eqs. (|2.2ip and (12.56)) 
are functionally identical, the same three forms of Kubo 
formula Eqs. (|3.2j|3.5p for the response of T to Hi hold 
even in the presence of a magnetic field, provided one uses 
the appropriate strain generator as given in Eq. (|2.49[) . 

The different forms of the Kubo formula can be viewed 
as related by use of Ward identities. Generally, Ward 
identities are the consequences of symmetries or conser- 
vation laws for response or correlation functions. In our 
case, the key relation, Eq. (|2.2ip or (|2.56p . was obtained 
by expanding the Fourier-transformed continuity equa- 
tion to first order in wavevector. 

There are some technical points about the derivation 
to discuss. These center around the assumption of time- 
translation invariance, which is usually assumed to hold 
for response functions such as eq. (|3.2p , on the basis 
that the unperturbed ground state is an eigenstate of 
the unperturbed Hamiltonian (or similarly at non-zero 
temperature, because the thermal weighting is station- 
ary). As we have mentioned already, for the systems we 
consider, in zero magnetic field a normalizable ground 
state is generally not available unless a confining poten- 
tial is included; this is because the system is in infinite 
volume, so that the strain generators can be defined. In 
Appendix [Cj we show that the preceding Kubo formulas 
are certainly valid as written in the presence of non-zero 
magnetic field, using the formalism of Sec. 111.21 and also 
more generally using the formalism of Appendix [AJ 

For rotationally-invariant systems, in which Hq com- 
mutes with angular momentum (the antisymmetric part 
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of Jfiu), the symmetries of X can be read off from the 
formulas. First, T^ v is symmetric, and therefore X^ va p 
is symmetric under /i <-> v. Next, in the stress-strain 
and stress-stress forms the part of X antisymmetric un- 
der a •<-> /3 also vanishes, using the assumption that the 
ground state is an Hq eigenstate (more detailed or care- 
ful arguments can be given along lines discussed at the 
end of Sec. IIII.51 and in Appendix [Cl below) . This holds 
without assuming the ground state is an angular mo- 
mentum eigenstate. Finally, for the strain-strain form, 
at first sight it may be less obvious that the part of X 
that is antisymmetric under fi v vanishes. However, 
for these components one can see that the time-integral 
term cancels the contact term (and similarly for the part 
antisymmetric under a -H- /3), again without assuming 
the ground state is an angular momentum eigenstate. 



Intensive form of the response and periodic 
boundary conditions 



III.2. 



The stress-stress form, eq. (|3.4[) . of the response 
X^ a p{ui), which is an extensive quantity, can also be 
understood in another way, as the response to a change 
in A using the Hamiltonian H\ in the x variables. In 
the rotationally-invariant case, this is the same as the re- 
sponse to a change in metric, and this point of view may 
be familiar to some readers. Using this approach, we can 
also obtain a Kubo formula in terms of intensive quan- 
tities only, which is compatible with periodic boundary 
conditions, and will be useful later. 

To rederive eq. (|3.4p . we work in x variables, and recall 

fx) 

that the integrated stress tensor we are using is 7) t „ , eq. 

(|2.31[) . and its expansion to order A in the x variables was 

given in eq. (|2.32l) (and the same forms are valid with a 

fx) 

magnetic field). Then the linear response of T^' to the 
strain, divided by —iuj + so that this is actually the re- 
sponse to dXap/dt, and finally dropping the distinction 
between x and X variables (as we require only the linear 
response) , gives exactly eq. (|3.4p . The contact term came 
from the expansion of the stress to order A in eq. (|2.32p , 
just like the familiar diamagnetic term in the conductiv- 
ity response comes from expanding the current operator 
to order (A^ being the perturbing vector potential, 
the response to which is conductivity). 

A similar approach works for the linear response of the 
local stress tensor r^^(X) = (det A) _1 A QA1 A^r Q( g(x) to 
the uniform strain [the (detA) -1 factor is present be- 
cause this transforms as a density, or formally because 
the (5-function should be written in X space rather than x 
space]; the result is the same as eq. (|3.4p . but with 77^(0) 
in place of T^ u . Here we have used the same boundary 
conditions as for the previous derivation, that is, a sys- 
tem in an infinite volume. But a similar derivation also 
works for periodic boundary conditions. These are de- 
fined in x variables as periodic boundary conditions on 
a box (or "unit cell") of fixed shape and size (indepen- 



dent of A), say a cube. The Hamiltonian is H\ as before, 
except for the different boundary conditions, and unin- 
teresting changes to the interaction potential to ensure 
it is periodic. When a magnetic field is present (in two 
dimensions), the flux through the unit cell is fixed in- 
dependent of A also, and must be an integer number of 
flux quanta (the flux quantum is 2ir in our units); use 
of the operator B is not required here. In X variables, 
the box has periodic boundary conditions described by 
A, as X = A T x, while the Hamiltonian is indepen- 

dent of A, except possibly through the interaction po- 
tential, the boundary conditions, and also the magnetic 
field strength varies with A so that the flux through the 
unit cell stays fixed. We note that with these boundary 
conditions, translational invariance holds strictly. Eq. 
(I2.32p does not hold due to the non-existence of J^ v in 

a finite-size system with these boundary conditions, but 

fx) 

there is still an expansion of t\ w (X) to order A in the x 
variables: 



r (X) 



(X) 



r (0) 

///,' 



(x) 



d\ a f 



\ a p + 0{\ 2 ). (3. 



A=0 



Another benefit of these boundary conditions is that 
normalizable energy eigenstates always exist. Then the 
Kubo formalism, working in x variables, leads to the re- 
sult for linear response of the stress at x = in the 
ground state (or at nonzero temperature) to a uniform 
rate of strain, which we call 



ioj 




dX 



a/3 



A=0/ 



H J dt I d d xe 



^ ) (O,t),rS(x,0) 



3.7) 



the integral over x is restricted to the box. (There is of 
course also a similar formula for X with these boundary 
conditions; one expects the leading, extensive part of X 
to be independent of the choice of boundary conditions.) 
Now we can take the infinite size limit (with particle den- 

sity held fixed, as always). In the limit, because t^ v '(0) 
is a local operator, its expansion to order A Q ^ in terms 
of the commutator with J a p is again valid. Then finally 
we have [dropping the superscript (0)] 



+ j dt J d rf a; e ^ + *([T^(O,i),T Q/3 (x,0)]) o |.(3.8) 

This agrees with the argument sketched just above, which 
began with TV particles in infinite space. 

We can also transform this stress-stress expression into 
a stress-strain form. As we have already passed to the 
infinite system, we can use the relation eq. (|2.2ip and 
integration by parts again, to obtain 

/>00 

X^ap(u) = -i dte iu+t {[T„ v {Q,t), J a p{0)]) . (3.9) 
Jo 
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Note that a strain-strain form for the response is not 
available in terms of intensive quantities. 



III. 3. Viscosity from the response function 

It is natural to ask how the extensive and intensive 

forms are related. Xnuap(u) is not simply X^ va fi{iS)l L d , 

because of the following simple point. Recall that X is 
fx} 

the response of (Tp, v ) to strain, while x is the response of 
fx) 

( T vu> (0)) instead. Assuming a homogeneous fluid state, 
these should be related by 



(^> = ^(i x) (o)) 



(3.10) 



where L d means the volume of the system, which is the 
A-dcpcndent volume of the box. L = L |A=jdetA, cal- 
culated either in the X variables as the volume of the 
A-dependent box using the fixed, standard Euclidean 
metric, or in x variables for a fixed box but with the 
A-dependent metric. (Actually, without rotation invari- 
ance, these metrics might not be uniquely defined by our 
models, however the formula is still correct, because the 
important point is how the volume varies with A, which 
is always through det A only.) Then for the response at 
first order, we obtain 

X 'nvapiw) I 'L d = Xfiisa0(u) ^S^SapP, (3.11) 

where in the last term we use the result at zeroth or- 
der in the strain, that in the thermodynamic limit the 
expectation of the local stress is 



<r$(o)> = <WP, 



and this becomes 5 M „Pi n t in the presence of a magnetic 
field in two dimensions. 

The response function x 1S almost exactly what we 
need to obtain the viscosity, which is supposed to be the 
local stress response to a uniform dX^/dt, and so is an 
intensive quantity. However, the expectation value of the 
stress that we just discussed will respond to static (time- 
independent) strains — see for example the hydrodynamic 
forms (response local in space and time) in eq. (|1.2|) . As 
X is the response to the time-derivative of the strain, 
these elastic moduli at zero frequency will show up as 
the coefficients of singularities ~ i/uj + in x- If non-zero, 
they might be confused with viscosity coefficients that 
happen to diverge at u = 0. 

To deal with this and obtain expressions for the vis- 
cosity tensor at all frequencies, we will remove from x 
these static, and thus equilibrium, elastic coefficients, by 
subtracting the zero-frequency value of the response to 
strain, that is the elastic moduli divided by iuj + , and 
call the remainder the viscosity tensor. For the homo- 
geneous fluids we consider, the expectation of the stress 
obeys eq. p.!2[) . which has an obvious generalization in 
the presence of a static strain, and is affected only by a 



dilation. Then we arrive at our definition for the viscosity 
tensor at frequency u in zero magnetic field, 



i ( dP \ 



(3.13) 



or, in terms of the inverse compressibility k 1 
-L d (dP/d(L*)) N , 



-Sn v S. 



a/3- 



(3.14) 



This, along with Eqs. (|3.3H3.5|) and ([3~TT]) . or Eqs. (pTgjt 
and (|3.9p . give the Kubo formulas for viscosity in the 
absence of magnetic field. 

The case of two dimensions with a magnetic field is 
similar, but the dilations involve a rescaling of magnetic 
field, as we have seen. Then we arrive at 

T)ituap(u) = X^vafj{^) + 5 ^ 5 a p L 2 ( ^(£2) ) 



IK; 



CO' 



(3.15) 



where n ; . 



-L 2 (dP int /d(L 2 j) v N is the "inverse in- 
ternal compressibility" , with the second partial deriva- 
tive taken with BL 2 (or the filling factor) held fixed. 
If we write the energy of the system as E(N,L 2 ,B) — 
L 2 e(y, B), where e(v, B) is the energy density as a func- 
tion of Landau level filling factor v — nfio/B, where 



N/L 2 is the particle density and 



(3.12) flux quantum 



27r in our units), then 



hc/e is the 



P iat (u,B)=B[ d£ ^' B) ) -s(u,B), (3.16) 



dB 

32, 



K -Hu R\ R 2 ( d e ^ B ) 



(3.17) 



In the fractional quantum Hall effect, one encounters 
incompressible fluids, which have vanishing compressibil- 
ity at zero temperature. This refers to the usual com- 
pressibility, which can be related to the change in den- 
sity with chemical potential with the magnetic field fixed, 
and so differs from the internal compressibility consid- 
ered here, which is well-defined and usually non-zero and 
finite. (Later, we will see that the latter can be extracted 
from the q 2 part of the conductivity.) 

Further, there are similar results at non-zero tempera- 
ture. For the (internal) compressibility, there is the ques- 
tion of what is held fixed in taking the partial derivative 
of the pressure: the temperature or the entropy. We will 
discuss this point in Sec. IIII.41 an d argue that it is the 
entropy, as mentioned in Ref. Il4l After this point, our 
discussion applies to both zero and non-zero temperature, 
unless otherwise noted. 
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III. 4. Non-divergence of bulk viscosity 

Next we comment on the subtraction of terms contain- 
ing the thermodynamic inverse (internal) compressibility 
in the expressions for the viscosity in terms of the re- 
sponse function, Eqs. f|3 . 14[) and (|3.15p . in either the ex- 
tensive or intensive forms. We will be led to the striking 
conclusion that the bulk viscosity is finite, not infinite. 
Consider first the case of zero magnetic field in d dimen- 
sions, and the extensive response X in finite size. We 
can show in general that the leading contribution to the 
diagonal or trace part = X w „„ / d 2 of the response 
function as lo — > is given by the derivative of PL d with 
A. In the stress-stress form, X<\ is 

X d (u) = ^ T {[T(0),K(0)}) 

dte^ +t ([T(t),T(0)}) Q , (3.18) 

where we have defined T = T^/d with d the dimension 
of space. Asymptotically as u — > 0, 

1 



u + X d ~~{[T(0),K(0)}} 



lim / dt 



[T(t),T(0)]) . (3.19) 



This is the response of the system to a static dilation at 
zero wavevector. By inserting a complete set of energy 
eigenstates |e) of H into the commutator in the second 
line, we obtain 



1 



[T(0),K(0)} 



2i 



I Tie) 



En — E P 



(3.20) 

(the restricted sum over states is really a principal part, 
arising from careful use of e — > 0). As explained already 
in Sec. IIII.21 working in the x variables the first term on 
right-hand side is the expectation of the change in the 
operator T( x ) under a dilation. Then time-independent 
perturbation theory using the x variables shows that the 
right-hand side is the change in the expectation of 
due to the dilation, the second term being due to the 
change of the ground state. Thus 



Lo+X d 



\d(L*)) N 



. ( d{PL d ) 
tin' 1 -P)L d , 



N 



(3.21) 



showing that the response to a dilation at leading order as 
lo — > is given by the difference between the inverse com- 
pressibility and the pressure — it is an elastic response. 
This is precisely the term we subtract from X to get the 
viscosity in Eq. (|3.13|) . We conclude that there can be no 
divergent bulk viscosity at zero frequency. The argument 
goes through similarly in the presence of a magnetic field 



in d — 2, with P and k _1 replaced by P mt and re- 
spectively. It also goes through similarly for the intensive 
formulation of viscosity in terms of \, eq. (|3.7[) . 

This argument may give rise to some unease. The shear 
modulus of a fluid is of course zero, while the preced- 
ing argument may lead us to expect that it is given by 
a similar time-independent perturbation theory expres- 
sion, which (as we will see) is non-zero in some cases. 
This must be related to the fact that the conventional 
way to obtain the elastic moduli (including the inverse 
compressibility or bulk modulus), which are susceptibil- 
ities, from the stress-stress response function, would be 
by taking the thermodynamic limit of the (intensive) re- 
sponse first, then taking to — > before q — > 0. But we 
have q = 0, then took lo — > 0, and finally the thermody- 
namic limit. 

One can find susceptibilities directly at q — 0, pro- 
vided one is careful with the order in which one takes the 
derivative with respect to strain and the thermodynamic 
limit. The usual "thermodynamic" formulas involve fi- 
nite size, but the size is treated as large and discreteness 
effects are ignored when differentiating the ground state 
energy, which corresponds to taking the limit first. (Here 
and in the remainder of this section we concentrate on 
the case of zero temperature.) The perturbation theory 
formula for the derivative (in finite size) , as in the right- 
hand side of eq. (I3.20[) . takes the derivative of the ex- 
pectation in a state that varies continuously with strain. 
If the ground state energy level does not cross others as 
the strain is varied, the result is most likely independent 
of the order of the limit and the derivative. We expect 
this is the case for pure dilations, leading to the pressure 
and inverse compressibility. But perturbation theory can 
break down if energy levels cross, especially if they do so 
on a set in A space that becomes dense as the size goes to 
infinity. Then if by "ground state" we mean the lowest 
energy state for given strain (as in the thermodynamic 
formulas) , this state vector changes discontinuously with 
strain, and this effect cannot be picked up in perturba- 
tion theory. Taking the derivative after the limit will 
give a different result from taking it before, at whatever 
strain it is taken. This is the case for shear strains, in 
some gapless systems. 

As an example, consider the free Fermi gas. Using pe- 
riodic boundary conditions on a rhomboid-shaped box, 
we can examine the expectation of Hq and of T^, and 
their variation with shape to first order. In X variables, A 
describes the shape of the box, and the metric in Hq and 
in T^ v is the standard one. Then as is well-known, the 
many-particle ground state for a given box is constructed 
by occupying all single-particle plane-wave states with 
wavevectors k inside the Fermi sphere; the radius of the 
sphere is chosen to obtain the correct particle number TV. 
Under shear and dilation, the k points move around, and 
every so often points enter or leave the Fermi sphere. 
Under a pure dilation, the volume of the Fermi sphere 
changes so that exactly the "same" set of k points (up to 
a rescaling) is always occupied, and no levels cross. But 
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during a nonzero shear (or shear and dilation together), 
k points do enter or leave the Fermi sphere, and so the 
ground state energy levels cross. In the limit, the ground 
state energy density depends on the particle density in 
the system, but not on its shape. Hence derivatives of 
this energy density with respect to shear vanish, as ex- 
pected. Derivatives with respect to dilations give the 
pressure and inverse compressibility, and because no lev- 
els cross as the scale factor det A is varied, the same result 
would be obtained by differentiating the energy of the fi- 
nite volume system, and taking the limit afterwards. But 
for derivatives with respect to shear, in finite size the sec- 
ond derivative is nonzero, though it becomes undefined 
on a set of measure zero at which ground state levels 
cross. Through the above formulas, this leads to an in- 
finite shear viscosity, while the bulk viscosity is zero (we 
give further details for the free Fermi gas in Section IV. II 
below). We expect, though we do not have a rigorous 
proof, similar behavior for an interacting Fermi gas in a 
Fermi-liquid phase, so that the pressure and compress- 
ibility can be obtained by using dilations either before 
or after the thermodynamic limit. This means the i/uj + 
terms in X^/L d or the corresponding part of x are can- 
celed by the subtraction, and the bulk viscosity cannot 
be infinite. But the shear viscosity will be infinite at 
zero temperature, in agreement with calculations of its 
temperature dependence^. 

For the first derivative of ground state energy with re- 
spect to shear, that is for the traceless part of the expec- 
tation of the stress, the situation is slightly different from 
that for the second derivative. Taking the limit first, the 
ground state energy density will be independent of shear. 
If, as an interval of a path in A space is traversed (with 
tr A fixed), there are no level crossings, then the limit of 
the first derivative along the path will also be zero. But 
suppose there are level crossings on this path, and the 
spacing of them goes to zero in the limit, so that the 
set of positions of level crossings is dense (as for the free 
Fermi gas). As the derivative taken after the limit is zero, 
the lowest energy levels as a function of position on the 
path must be close to a sequence of overlapping parabo- 
las, with the minimum of each at the same energy. We 
have seen that the second derivative is of order one, and 
so the first derivative on any of these curves, anywhere 
within the interval in which it is the lowest, will tend to 
zero in the limit, because the distance in A from the min- 
imum of that curve goes to zero, as the level crossings 
become dense. So (except on a set of measure zero in A 
on which the first derivative is not defined) the traceless 
local (or intensive) stress in a fluid state does go to zero 
in the thermodynamic limit, as claimed earlier. 

At non-zero temperature, one can make different but 
related arguments. In particular, the trace part uj + Xd 
is now identified as — i(P — Kg 1 )^, or similarly with a 
magnetic field in two dimensions, where the inverse isen- 
tropic compressibility, 1 = -L d (dP/d(L d )) Ni s 

at fixed entropy S (and also fixed v for the isentropic in- 
ternal compressibility). This is because, as for the pres- 



sure, the response function can be identified as the par- 
tial derivative taken under the thermal average, with the 
probabilities and hence the entropy, held fixed. These 
compressibilities must be non-negative for stability. The 
appearance of the isentropic compressibility in this limit 
of the response is frequently obtained from hydrodynamic 
considerations, rather than directly from the stress-stress 
response, as here. 

In the preceding arguments, we have taken the fre- 
quency to zero in the response function before the ther- 
modynamic limit. However, we will be using the intensive 
functions X/L d , x and r\ in the thermodynamic limit, and 
the behavior of these as the frequency tends to zero sub- 
sequently. The different order of limits does not appear 
to be a problem. For example, in the zero-temperature 
case, we can use the intensive form in x variables to study 
the response of the trace of r to a low-frequency dilation. 
The leading part comes from the ground state adiabat- 
ically following the dilation, and gives the inverse com- 
pressibility. The real part of the bulk viscosity involves 
transitions to excited states of the unperturbed system, 
for which the available phase space is usually small (as 
in a Fermi liquid, for example) or zero. Hence we expect 
that at zero temperature, in general the real part of the 
bulk viscosity actually goes to zero at zero frequency. 

III. 5. Spectral density, sum rules, and positivity 

We can derive a spectral density for the viscosity ten- 
sor, and a sum rule for it, by following a standard 
method, starting from the convenient stress-strain for- 
mulation, in either extensive or intensive forms. For the 
extensive form, using Eqs. (|3 . 3[) . we define the spectral 
density function byi£ 

/oo 
dte^([T^(t),J a0 (O)}). (3.22) 
-oo 

In many cases of linear response theory, such a func- 
tion would be the imaginary part of the corresponding 
retarded response function, although for transport func- 
tions, such as conductivity as well as viscosity, the divi- 
sion by — iu (in the current-current, respectively stress- 
stress, forms) means that the spectral density is propor- 
tional to the real part of the conductivity, if certain sym- 
metries such as time-reversal and reflection symmetry are 
unbroken. But this is not generally the case: X" is not in 
general reat^, but does consist (in a finite size system) of 
a sum of (^-functions in w, with tensor- valued coefficients. 
There is a spectral representation, 

/oo V" (in') 

du> 7* M . (3.23) 
-oo - U)' 

This shows that if X" is real, it is the real part of A, as 
expected. X^ va ^{ui) is analytic in the upper-half com- 
plex oj plane, and there are also corresponding Kramers- 
Kronig relations between X" and the complementary 
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P art X' = [Xp Ua p(u) - x 'J va/3 (u))]/i (for real u), 

which is the imaginary part of X when X" is real. Fi- 
nally, the definition leads immediately to a sum rule for 
the total spectral density, 

— X'^ aP { U ) = -i{[T^(Q), J a p(0)]). (3.24) 

The "sum" on the right-hand side is real, because the 
expectation of a commutator of self-adjoint operators is 
imaginary. The sum rule can also be viewed as describ- 
ing the lu — > oo limit of X^ a p(uj), using the spectral 
representation on the one hand, and integration by parts 
from the stress-strain definition of X, together with the 
Riemann-Lebesgue lemma, on the other; this explains its 
relation with the contact term in the stress-stress form. 
There are related results for Xnvap and for (Sim- 
ilar sum rules were also discussed in Ref. [Til but only in 
the absence of Hall viscosity.) One would expect that the 
right-hand side is symmetric under \xv f-> a/3, which can 
be shown in certain limits, as we discuss below. That is, 
the Hall viscosity cancels from the sum rule. This would 
then be similar to the case of conductivity, in which the 
Hall conductivity cancels from the sum rule. 

For further arguments, the stress-stress form of the 
spectral density is most convenient. The only complica- 
tion here is the factor 1/lo + in the formula for X, which 
causes the appearance of 5(u) terms in the spectral den- 
sity, in addition to (5-functions that come from the stress- 
stress time-integral term. The former correspond to the 
terms discussed in Sec. 1111.41 We can obtain the spec- 
tral density by multiplying X" in eq. (I3.22[) by iu), which 
means differentiating with respect to t under the integral, 
and then using relation (|2.2ip once again (after shifting 
the t-dependence onto J a p). This has the effect of remov- 
ing any 5(u) terms from X" , one of which we know is the 
inverse compressibility term. Reinstating these terms, 
one has 

i r°° 

+ 2^J tee^(\T^(t),T a pQS)]) (3.25) 

(in finite size, the time-integral expression should not 
produce any 5-function at oj = 0), where the constant 
tensor C^ap is 

C/wap = - i ([7)^(0), J a p(Q)]} 

/>oo 

-i ^e- £ *([T^(t),T Q/3 (0)]) o , (3.26) 
Jo 

and is real. Then, using arguments presented in e.g. Ref. 
flEl . one can almost conclude that X"(uj), viewed as a 
matrix with rows and columns indexed by the pairs \xv 
and a/3, respectively, should be Hermitian, and also pos- 
itive semidefinite, for all u>, and that at zero frequency 
it (i.e. C^ va p) should be real. We say "almost" because 
Forstcr's discussion 15 does not include the contact terms 



in our stress-stress form, which contribute to X" only 
at zero frequency. Hermiticity holds at non-zero fre- 
quencies, and the time-integral term in C^ ua p is real 
and symmetric, so to obtain symmetry of the matrix at 
zero frequency, we would need the contact term coeffi- 
cient — i{[T^ u , J a @]) to be symmetric under the exchange 
[iv O a/3. This is not yet obvious in general, and we dis- 
cuss it further below. For positivity, the energy absorp- 
tion argument^ does allow for contact terms, and those 
that occur in the positivity statement at zero frequency 
are automatically symmetrized, as they arise from the 
second derivative of the Hamiltonian with respect to the 
perturbing field, in our case of H\ with respect to A. The 
part of C^ va p8{ui) that vanishes on taking the trace on 
fiu and on a/3 would either be zero, or would represent 
an infinite shear viscosity, as discussed in Sec. IIII.4I At 
zero frequency, the infinite shear viscosity (if any) must 
be positive. For arbitrary frequencies, in simple cases 
such as with time-reversal and reflection symmetries, in 
which case the X" matrix is symmetric, these conditions 
imply that the real parts of the shear and bulk viscosi- 
ties must be positive (like the real part of the symmetric 
conductivity tensor). In general, at non-zero frequency, 
the condition that the Hermitian matrix X" be positive 
semidefinite involves the imaginary part of the Hall vis- 
cosity, and not only the shear and bulk viscosities (just 
like the case of conductivity, in which the imaginary part 
of the Hall conductivity enters) ; note that the discussion 
in Sec. [J was for frequency- independent or zero- frequency 
viscosity coefficients only. At zero frequency, one state- 
ment of positivity is for the intensive function x"( w )j 
which differs from X" by the pressure term in the trace 
part; it implies that the inverse (internal) compressibility 
should be non-negative, in agreement with a consequence 
of thermodynamic stability. On the other hand, positiv- 
ity of the trace part of X" implies that — P > 
as well. This can also be viewed as a consequence of 
thermodynamic stability, by using tr A as a "generalized 
coordinate", rather than the volume L d . 

If the zero-temperature real part of the shear viscos- 
ity is a (5-function (that results in part from the contact 
term), then at non-zero temperature in an interacting 
system in the thermodynamic limit, the shear viscosity 
at zero frequency will generally be neither zero nor infin- 
ity, and as a function of frequency the (5-function becomes 
broadened. In this case the (5-function is cancelled by a 
contribution from the time-integral term, which also pro- 
duces the broadened peak (which tends to a 5-function 
at zero temperature). For such cases, the spectral rep- 
resentation above for finite size should be rearranged so 
there is no (5-function in the traceless part. 

The symmetry of the contact-term coefficient 

-iiW^Jap]) (3.27) 

under [iv O a/3 can be shown under some conditions or 
in some limits. Using relation (|2.21l) again, the antisym- 
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metric part is — | times the expectation value of 

[[Ho, JiMu],J a p] — [[Ho, J a @], Jpv] = [Ho, [Jfiu, Ja@]] 

(3.28) 

(this also gives the difference between the contact terms 
that would result in the x variables from use of T pv ver- 

sus T^ v at the end of Sec. III.1[) . The commutator of 
the two J's always contains only traceless parts, as the 
commutation relations can be rewritten 

— dfift^Jav ~^OLuJ n i r y) 3va(Jfj,fi ^^^^77/* 

(3.29) 

Taking the expectation value, and using the relation 
(|2.21j) . we arrive at two expectations of traceless parts of 
T pu , which we have argued are small compared with L d in 
the thermodynamic limit (or zero, in the magnetic field 
case in two dimensions). (In the rotationally- invariant 
case, we look at parts symmetric under fi v and under 
a -s-> /3 only, and then the commutator gives the angu- 
lar momentum, and the expectation is zero even in finite 
size.) Hence, if we are interested in the viscosity of the in- 
finite (thermodynamic limit) system, then the symmetry 
under [iv ■<-> a/3 does hold, and has the pleasing conse- 
quences mentioned above. 

One might wish for more in the finite-size case. Our 
basic definitions can, and perhaps should, be modified to 
make the contact term in the stress-stress form symmet- 
ric, by taking its symmetric part. It is not clear to us 
fundamentally why that would have to be done, but it 
would be useful anyway when performing numerical cal- 
culations at finite size, to remove the antisymmetric part 



that should not be present in the limit. Alternatively, 
perhaps there is some physical meaning to the antisym- 
metric part of the coefficient of S(lu) in X", even though 
it does not contribute to energy absorption, or to the 
physical (infinite-size) viscosity. 

IV. RELATION BETWEEN VISCOSITY AND 
CONDUCTIVITY 

In a Galilean-invariant system with particles that all 
have the same charge (which is one in our units) and 
mass m, the momentum density is m times the number 
current density. We will now use the intensive form of the 
stress-strain response function to derive a general relation 
between viscosity and conductivity for this case. We note 
that this does not require rotational invariance, provided 
that the inverse mass tensor is the same for each particle. 
In this section we will usually treat the zero and nonzero 
magnetic field case together, and so specialize to d = 2 
dimensions; other cases are handled similarly. 

We begin by using the translation-invariant system in 
a box of volume (area) L 2 with periodic boundary con- 
ditions, with Hamiltonian Ho- The continuity equation 
for momentum density, eq. (|2.51|) . can be written as 

^ Su\ + u c e uX ^J 9x(x,t) = d p T p „(x,t), (4.1) 

where lj c = B/m is the cyclotron frequency. (For 
B > 0, the following derivation still holds if there is an 
anisotropic mass tensor, by making a suitable modifica- 
tion of Lo c e^\ to a different tensor.) Due to the uniformity 
of the magnetic field, we can use this to derive additional 
modified Ward identities, in particular starting from the 
stress-stress retarded response function: 



q x q p I die 

10 



d 2 xe-^([T X „(x,t),T p0 (O,O)}) o = 



{iu5 vX +u c e vX ){-iu)5p p + oj c e ?p ) / dte lu+t I d'xe' 1 ^ ([g x (x,t),g p (0,0)]) 



tq^J d 2 xe-^ x ([r 7 ,(x), ff/3 (0)]) 
{iu)8p p -uj c ep p ) / Ae- iq,x ([g„(x),^(0)]) . 



(4.2) 



Expressions similar to this have been used for zero mag- right-hand side of Eq. (|4.2[) is the same as the one appear- 

netic field (w c = 0) by many authors, in particular Taylor ing in the Kubo formula for the electrical conductivity 
and RanderiaM, however even in that case our approach 
differs in some of the details. 

Assuming Galilean invariance, the number current 
density j(x, t) is related to the momentum density by 
j(x, t) = g(x, t)/m. Hence, the retarded function on the 
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(current-current response function), 



0V/3 (q,w) 



mui 

oo 



— / die 



'-' ' 1 d 2 xe^ x ([>(x,0,jMO,0)]) . 

(4.3) 



Meanwhile, the left-hand side of Eq. (|4.2I) includes the 
time-integral part of the intensive form of the response 
function for viscosity, with two factors of q contracted 
into it, presently in finite size, as in eq. (|3.7[) . We will 
now account for the remaining terms on the right-hand 
side, and aim to take the thermodynamic limit and only 
then expand in powers of q to second order, to obtain 
the response function Xnvap{w) from the conductivity. 
It follows from the preceding derivation that, as the left- 
hand side of Eq. (|4. 2[) is of order q 2 , all terms of order 
one and of order q on the right-hand side must cancel. 

First, still in finite size, in the last term on the right- 
hand side of Eq. (|4.2[) , we can use translation invariance 
to introduce integration over a variable x', divide by L 2 , 
and then evaluate, giving 

-L [ d 2 xd 2 x'e-^( x - x ')([ 5 ,(x), 5p (x')]) 



EK[{<' e " q ' x< }'H eiq ' x< }])o 



1 NT- 1 

L 2 

1 



= J2 { iNBe "P - 1" ( G p)q - IP ( G ^)o) > 

where G = 7r ! is the total momentum. The ther- 
modynamic limit of this exists, and contains no term of 



order higher than q as q —> 0. 

Next, the second term on the right-hand side of 
Eq. (|4.2|) can be manipulated to produce the contact term 
in x- First, we take the thermodynamic limit, and then 
the part quadratic in q, by expanding the exponential 
e -iq-x £ £ rs j. orc j er w e a i so use translation invariance 

again, so that the expression contains [t^^O), gp{— x)]. 
Then we recognize the occurrence of the first term of the 
strain generator, if we write Eq. (|2.49[) in the form 

Jjj,v J d X X ^Cji/ (x) -)- ^ g ^€ ua X ^X a 



+ -5^{l3,E({x%13)} 



(4.4) 



(or similarly without the terms containing £>, if B = 0). 
Then we can cast the contact term from Eq. (|3.8p in the 
form, 

1 1 f 

—r <[t m „(0), J a p]) = / d 2 x ([T„ v (O),x a g (x)]) Q , 

OJ OJ J 

(4.5) 

where for B ^ the last two terms in Eq. (14.41) do not 
contribute, because they fail to commute only with the 
kinetic part of the stress tensor, and the result of that 
commutator always contains a product of a particle co- 
ordinate with a delta function of that coordinate at the 
origin. This means that the term in question can be ab- 
sorbed into the time-integral term on the left to produce 
precisely x contracted with two qs. 

Hence, Eq. (|4.2p turns into a relation between the 
conductivity and a symmetrized part of the intensive 
strain-stress response function; we define Xpuapi 10 ) — 
\\Xnvap{u) + Xavnp{u)\, then 



<9 2 er A p(q, uj) 



dqudq a 



(4.6) 



q=0 



This is the central result of this section. When the magnetic field is zero, one can simply substitute uj c = in the 
above equations. Let us note that a similar relation holds in three dimensions: 



<9 2 CTA P (q,u;) 



dqp,dq a 



(4.7) 



q=0 



where b = B/|B| is a unit vector in the direction of the 
magnetic field. 

These results lend themselves to a simple interpreta- 
tion. For simplicity, we will concentrate on two dimen- 
sions. If we expand the wavevector-dependent conduc- 
tivity in powers of the wavevector q, 



cr Ap (q,cj) = o- A pH + oipfa 



(4.8) 



electric field) is fixed by Galilean invariance to be 



TO (u>+ 2 — UJ 2 ) 



(iuj + S Xp - w c e Ap ) (4.9) 



By Eqs. (|4.6[) and (|4.9[) . the second order in q term is 
related to the viscosity through 



then the zeroth order term (the response to a uniform cr Xp( c l^ L0 ) ~ a \v ( u ) —IpXp^afii^) qa&p p ( w )- (4-10) 
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One can intuitively understand this expression as follows: 
in the presence of a nonuniform electric field E(q, u>), to 
the leading order the system responds with a nonuniform 
current, (jp(q,u;)) = a^l (oj)E\(q, u)), which implies a 
strain rate u a p (q, w) = {iq a )<j^{uj)E\(<\, oj)/n (sym- 
metrization over a and /3 is not important, as Xnvapip) 



is symmetric with respect to these two indices). The 
strain rate results in an average stress (r M1/ (q, oj)) = 
~Xpvafi (w)u Q| a(q,w), whose spatial derivative gives a 
contribution to the effective electric field acting on the 
particles, Ef, s (q, oj) = —iq fJ ,(T IMl ,(q,u}))/n, which then af- 
fects the current (and the conductivity) via cr^(w), re- 
sulting in Eq. (l4"T0|) . 



By Eq. (|3.14l) or (I3.15[) . to arrive at the viscosity one should subtract the inverse (internal) compressibility contri- 
bution from x^a/jM- Defining rj pua0 {u) = \[r]^ af} {uj) + Jfe^M], we have 



<9 2 cr Ap (q,w) 



dq p dq 



(uj5 



pP 



IK- 



q=0 



2uj+ 



(8fj,v6 al 3 + 5 P p5 va ) . 



(4.11) 



For rotationally-invariant two-dimensional systems, there are only three independent (frequency-dependent) co- 
efficients of viscosity: the bulk viscosity £, shear viscosity 7/ sh , and Hall viscosity r] H [cf. Eqs. (|1.6I) - (|1.7I) ]. The 
relation (|4.6I) can be used to extract the viscosity coefficients at all frequencies from the conductivity [taking q in the 
x- (i.e. 1-) direction without loss of generality], 



CM 



W" 1 



X1111M " X1212M {(^ ~ w c) kii(q,w) -CT 2 2(q,w)]} 



q=0 



V sh (u) = X1212H =~Y {w 2 CT 22 (q,w) + w 2 ern(q,w) + 2iujuj c a H (q,uj)} 



q=0 



Xm 2 M - XmiM _ m2 92 



IUJUJ C 



[on (q,w) +cr 22 (q,w)]} 



(4.12) 
(4.13) 
(4.14) 



q=0 



where cr ff (q,cj) = [cr 12 (q,w) — cr 21 (q,w)]/2 is the Hall conductivity. Here we have separated Q(uj) from wrj/w, 
according to our analysis of the trace part of the viscosity tensor in Sec. IIII.4I For zero magnetic field (lu c = 0, and 
k^J = kT 1 ), relations similar to the first two are fairly well known, but often are written in terms of the transverse 
and longitudinal parts of <j\ p ] see for example Ref. [lj, in which however the k _1 term is absent. The Hall viscosity 
decouples in this case. 

Alternatively, again for a rotationally-invariant system in two dimensions, we can invert Eq. (|4.6[) to obtain the q 2 



(2) 

part of the conductivity tensor, cr^ p (q, w). For vanishing magnetic field we find: 



2,,,+2 



CM 



qxQp + V sh { u ')'l 2 S\p + r] H (u))q 2 e\ p 



(4.15) 



For zero magnetic field, we see that the Hall viscosity can be obtained just from the q 2 part of the Hall conductivity, 
at all frequencies. 

In the presence of a magnetic field, the expressions are more complicated, and for later use we retain only terms 
that are non-vanishing at zero frequency, assuming that £ and rj H do not diverge as ui — > 0, and that r] sh does not 
diverge more rapidly than 1/lu + . Then wc have as u — > 



4 2 >M 



1 



C(w = o) + 



ip£p 



„xq a e ap + V sh {Lo)q 2 6x p + ( r, H (to = 0) - + ^rf\u) ) q 2 e Xp 



(4.16) 

This can be used to obtain r\ H (to = 0) from the q 2 part of the conductivity at uj — > 0. If r] sh is non-diverging as ui — » 0, 
the antisymmetric part of the equation reduces to the relation found by Hoyos and Son 8 between the Hall viscosity, 
internal compressibility, and q 2 part of the Hall conductivity at zero frequency in the presence of a magnetic field, 
which they obtained for a gapped quantum Hall system at zero temperature (see eq. (|3.17l) for the equivalence of k^* 
with their expression). That relation is now seen to hold even when £(w = 0) and i] sh (u) = 0) are non-zero but finite, 
with a generalization for the case of diverging 77 . Our formulas above give the generalization to all frequencies as 
well. 

I 



EXAMPLES OF THE KUBO FORMULAS 
AND CONDUCTIVITY RELATION 



tional invariance, at zero temperature except where oth- 



In this section, we consider applications of the Kubo 
formulas to various simple model examples with rota- 
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erwise stated. 



V.l. The free Fermi gas 

As a first application of our Kubo formula, let us com- 
pute the viscosity of the free Fermi gas in d dimensions 
at zero temperature. The unperturbed Hamiltonian for 
this system is 



(5.1) 



For a system of N particles enclosed in a very large box, 
the ground state energy E has the form 

E Q = (H Q ) K{L d )-i. (5.2) 

From Eq. (|2.21|) . the stress tensor is given by 

= -i [H Q , J MV ] = -Y»t- (5-3) 

Because T^ v is independent of time, the time-integral 
term in the stress-stress form of the Kubo formula Q3.4p 
vanishes, and we easily find 



1 



([T^(0),J^(0)]> 

(5^/3 {T vix )q + S^a (T v p} Q ) 



2i 



Eq (Sp/jSva + d^aSvfj) , (5.4) 



where we have used rotational invariance to express 
(Tfj,u) in terms ot Eq: P = 2E /(dL d ). Lastly, we have 
that 



N 



7 En 



N 



d 2 L d u ' 
(5.5) 

and so from Eq. (I3.14[) . we see that the zero-temperature 
viscosity tensor of the free Fermi gas is given by 

rf~+J^3 + &na?>vP ~ "-^fnvSapj . (5.6) 

This is purely a shear viscosity — the bulk viscosity of 
the free Fermi gas is not only not infinite, as anticipated, 
but identically zero at zero temperature. Furthermore, 
the dissipative shear viscosity coefficient r/ sh , given by the 
real part of the scalar prefactor of rj^ ua p , is 



h 2irE 
71 (W) = ^ (W) 



(5.7) 



which is zero for all w/0, and infinite at u) = 0. Such 
delta-function divergences in response functions are char- 
acteristic of noninteracting systems, and indicate that in 
response to shear strains, the free Fermi gas accelerates 



without bound. In an interacting Fermi gas at non-zero 
temperature, the ^-function is broadened and the zero- 
frequency shear viscosity is finite. As the temperature 
tends to zero, it approaches a 5-function; the limit is 
continuous in the space of distributions. 

To conclude this subsection we note that one can derive 
the same result by first finding the q-dependent conduc- 
tivity through elementary Green function techniques, 



2i E, 



dm 2 uj +3 L d 



2 (<Wq 2 + 2q v q p ) 



+ 0(q 4 ), (5. 



and substituting this expression into the zero mag- 
netic field version of the viscosity-conductivity relation, 
Eq. (|4.11l) [or, equivalently, comparing it with Eq. (14.151) ]. 



V.2. The integer and fractional quantum Hall fluids 

Let us now compute the viscosity for non-interacting 
electrons in an external magnetic field. The Hamiltonian 
H is 



Hn 



1 i i 
> 7T„7T„ 



(5.9) 



We can take our unperturbed state to have the lowest v 
Landau levels occupied, in a region of area L 2 , and so 
the ground state energy is 



E = nsuj c L 



2racj)Q 



-B Z L Z 



(5.10) 



where as before s = v/2 is minus the average orbital spin 
per particle. However, this is unnecessarily restrictive 
here, and at any temperature and any average filling fac- 
tor v the average energy is again (-ffo)o = nsu> c L 2 . From 
Eq. (|2.56|) . we have that 



We will calculate X using both the strain-strain and 
stress-stress Kubo formulas. It will be convenient to work 
in the symmetric gauge, where the dilation generator K 
takes the simple form 



(5.12) 



We can diagonalize the Hamiltonian Hq with two sets of 
creation annihilation operators. Writing z 3 = x 3 + iy 3 
and z 3 = x 3 — iy 3 , these are (see e.g. Ref. @) 



b l = 



(tT* + ) , 



6* - J** 
V 2 



(5.13) 
(5.14) 
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satisfying 



= 5. 



ij > 



= [&\a J ] = o. 



(5.15) 
(5.16) 



In terms of these operators, the Hamiltonian takes the 
simple form 



H 



B 



(5.17) 



The stress tensor can be written as 



Tfj, v — S^uHq- 



B_ 

2m 



(5.18) 

where t x and r z are the standard x and z Pauli matrices. 
The shear generator J*~ from Eq. (12.38)) takes the form 



a + a 



E^< 



r - 



(5.19) 



(note we could have written = ir^), and the dilation 
generator K can be written 



(5.20) 



The shear generators agree with those presented by Read 
and Rezayi^. The last ingredient we need is the commu- 
tation relations between V and the creation and annihi- 
lation operators. A computation shows that these are 



1 ,.it 



(5.21) 
(5.22) 



Now, using the stress-stress Kubo formula (|3.4|) . we 
find rather directly that 



£•(1 



— - [«W + ((5 /i( g(5 zyQ - e^e^) 



(5.23) 



where in uj c = B/m, B is again the value at the center of 
the narrow wavepacket over B values. Finally, applying 
Eqs. flOD and (|3~T71) . we have 



-Pint = nsw c 



En 



K inl = 2nSUJ c = 2 — 



(5.24) 
(5.25) 



and hence 



E 



' L 2 (uj+ 2 -4u%) 
-2w c (<5^ ct e /i/ 3 — 8npe a v) 



(5.26) 



At non-zero temperature, the result is the same, with 
the ground state energy density Eq/L 2 replaced by 
(Hq)o/L 2 . We notice that the trace on [iv or on a/3 van- 
ishes at all frequencies — the bulk viscosity is identically 
zero. The corresponding spectral density consists solely 
of ^-functions at the frequencies ±2w c , which represent 
transitions in which the Landau level index changes by 
±2. This is because of the quadrupolar nature of the 
stress, and the fact that only the operators b l appear in 
Tfj, u . In the uj — > limit, the viscosity reduces to the Hall 
viscosity 



Vi*vap(u = 0) = \ns (c)„ Q e M £ - <^e Q „) , 



(5.27) 



in agreement with known result o 3 ' 4 ' 18 . It is not surprising 
that the zero-frequency bulk and shear viscosities vanish 
when the temperature is zero and there is a gap in the 
spectrum; that they do so in other cases as well is due to 
the non-interacting nature of the system. 

To make contact with the adiabatic calculation of the 
Hall viscosity^— we can also calculate the viscosity 
from the strain-strain formula (|3.5[) . Naturally, this must 
give the same result as Eq. (|5.26[) . however it is enlight- 
ening to see how this comes about. If we evaluate just 
the equal-time contact term, we find 

-i (IV, Jap]) = | ( ^ - £ («' V + §) o ) 



(5.28) 



The intra-Landau level (a^a) term is larger than O(N), 
but we know that it must cancel. On calculating the 
time-integral contribution, we find as uj — > 

poo 

uj+ dte iu+t {[J^{t),J afi (0)]) ~ 



(5.29) 



Here the i/uj + term arose because there is a part of the 
correlation function that is linear in t for large t. That 
is present because the trace K = J w has time derivative 
— i[K, Ho] = —T^p = —2Hq, which is time independent. 
Combining these, we see that the intra-Landau level con- 
tributions to X exactly cancel in the final result. The 
traceless part of this result was obtained previously by 
considering the transport of degenerate subspaces in the 
disk geometry in the infinite plane in the adiabatic trans- 
port formulation of viscosity^. 

We can also use the strain-strain form of the Kubo 
formula (|3.5[) to calculate the uj — > traceless viscosity 



24 



of a fractional quantum Hall fluid in the lowest Landau 
level. Although the presence of interactions makes ma- 
nipulating the stress tensor rather cumbersome, we can 
calculate the necessary matrix elements of the traceless 
strain generators (|5.19[) for certain trial states. Because 
of the non-standard time dependence of the dilation gen- 
erator K, we avoid computation of the diagonal response 
function. 

Let us consider a fractional quantum Hall system in 
the lowest Landau level. We assume that the interaction 
is rotationally invariant, and commutes with all the Lan- 
dau level raising and lowering operators b' l \ b l . Further, 
we assume it is one of the "special" Hamiltonians for 
which exact zero-interaction-energy ground (and edge, 
and quasihole) states lying in the lowest Landau level are 
known (for more details on these, see Ref.©. bmcc we 
are only concerned with the traceless strain generators, 
we are free to work with states with a fixed magnetic field 
B. We denote by |0) the unique minimum angular mo- 
mentum ground state of the fluid. Using the strain-strain 
form of the traceless response function X sh in Eq. (|3.5I) , 
along with the commutation relations Eq. (I2.9[) . we have 

*#apH = S va (0|J*» - Stf (0| j£\0) 

+ oj + / dte** +t <0|[j£(t),J$(0)]|0>. 
Jo 

(5.30) 

Using the fact that the system is in the lowest Landau 
level, we can evaluate the first term above to get 

i 

+ u+ dte*- +t <0| [J#(t), J$(0)]|0>. 
Jo 

(5.31) 

In the limit uj — > 0, we can evaluate the time- integral 
term. Note that in this limit, the only nonzero con- 
tributions come from elements of the expectation value 
which are independent of time (there is no other kind of 
non-oscillatory time-dependence for these traceless com- 
ponents). Let {\D) : D = 0, 1, . . .} be an orthonormal 
basis for the subspace degenerate with the ground state 
(including the ground state |0)), and {|e) : e = 1,2...} 
an orthonormal basis for the space of all eigenstates with 
energy larger than that of the ground state (thus we are 
assuming a discrete spectrum, as in a finite system) . Now 
we use the following fact: Given a system with a discrete 
spectrum, and denoting by Po — J2d \D) {D\ the projec- 
tion operator onto the lowest-energy subspace, we have 

/>oo 

lim uj+ / dt e lu+t (0| A(t)B(0) |0) = (0| A(0)P B(0) |0) 

(5.32) 

for any operators A and B. (In fact, the identity contin- 
ues to hold if |0) is replaced by any state in the degenerate 



subspace.) Using this in Xf^ a g(u> —> 0), we obtain ex- 
actly the expression that results from adiabatic transport 
of a degenerate subspace as in Ref. 

X s * a/3 (u) = 0) = -i(0|VP^|0) +i(0\J a pP^ v \0), 

(5.33) 

where P± = 1 — Po- Thus it is the expectation of the 
commutator of J's, as in the contact term, but in the 
intermediate-state sum, the states degenerate with the 
ground state are omitted. (The noninteracting system 
considered above is a particular case, in which all intra- 
Landau level effects cancel; this does not occur in frac- 
tional quantum Hall states in interacting systems.) 

In the present case, using Eq. (|5.19|) and noting that 
J2i a i^ 2 leaves the ground state in the degenerate sub- 
space while a i takes it out of that subspace^, we have 

(o\J$#>J§ |o> - <o| J%P Jf v |0) = 

l - (0| a 1 2 af |0) (S^e av - 6 av e^) . (5.34) 
Inserting Eq. (|5.34[) into the Kubo formula (|3 . 5|) , we find 

\ i ij I ( 

= \~sN {5 va e^ - S^peav) , (5.35) 

(at leading order in the thermodynamic limit) where the 
second line follows from the matrix elements computed 
in Ref. 51 Thus again we obtain for the viscosity 

f] ll uap{^ -> 0) = \'sn{8 ua e f _ L p - 8^e av ) . (5.36) 



V.3. The viscosity-conductivity relation for 
quantum Hall systems 

In Sec. IIV1 we obtained a general relation between the 
viscosity and the q 2 part of the conductivity for a system 
in a magnetic field. We have already commented there 
that as lj — > we obtain a relation 

ip 2 |i^( q ,. = o)| q=0 = !7*( W = )-^+H^VV) 

uq x lo c oj c 

(5.37) 

which is more general than that of Hoyos and Son£, but 
which reduces to theirs [using eq. p.l7[l ] when the shear 
viscosity is non-diverging at u> — > 0. 

It is interesting to compare the full frequency- 
dependent relation with an exact calculation for the in- 
teger quantum Hall state at filling factor v, based on 
the results of Chen et al^ for cr Ml/ (q, uj). Extracting the 
quadratic in q terms in the conductivity from their cal- 
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culations we find: 



9 2 crii(q, uj) 



<9 2 <x 2 2(q,w) 



2/2 



d 2 a H (q,uj) 



q=0 



q=0 



uj c 2tt V w+ 2 - Auj 2 lo+ 2 - UJ 2 



Au: 2 



(5.38) 

icj+ 27r V w+ 2 - 4u 2 w+ 2 - oj 2 

(5.39) 



dg 2 



= -2- 



z^ 2 



q=0 



2tt Vw+ 2 -4w 2 w+ 2 



(5.40) 



where I — 1/vB is the magnetic length. Substituting 
these expressions into Eqs. (|4.12[) - (l4.14j) we arrive at (us- 



ing again K in J = vnuj Cl s = v/2): 



ri sh ((jj) =ns 
rj H (uj) =ns 



Alu 2 - oj+ 2 ' 
2uj 2 . 



(5.41) 
(5.42) 

(5.43) 



4w 2 - oj+ 2 ' 
in full agreement with Eq. (|5.26l) . 



V.4. Complex f-wave paired superfluids in two 
dimensions 



Lastly, we shall consider a complex £-w&ve paired su- 
perfluid of fcrmions in two dimensions. The model mean- 
field Hamiltonian we shall consider takes the form 

H = J d 2 x^{x) (-^ V ' - V<(x) 

+ J d 2 xd 2 x' A(x-x')^(x)^(x r ) + h.c, 

(5.44) 

where the pairing function A transforms as an ^-wave 
under rotations. We note that we are now working with 
a number-non-conserving system, in which the chemical 
potential appears as a parameter. Previously we worked 
with systems at fixed number, however, generalizations 
to fixed chemical potential (the grand canoncial ensem- 
ble) , or even as here to systems in which particle number 
is not conserved, should be reasonably self-evident. The 
terms in the Hamiltonian that violate number conser- 
vation also violate angular momentum conservation; the 
system is not rotationally invariant. However, the opera- 
tor Q = Cv^Jpu — -n^-N is still a conserved quantity, where 
^v^J^v is the angular momentum operator, and N is the 
number operator. 

Using the continuity equation (|2.20p . we find that the 



second-quantized strain generator is given by 
Jfiu = - J d 2 xx fl g u (x) 



I t . dip(x) 0<0 f (x) 
d x Xf, lip 1 (x) — — V(x 



dx v 



dx v 



(5.45) 



The general relations such as eq. (|2.21l) between stress 
and strain generators still hold, as do the Ward identities 
for the response function, and their consequences. We 
should view the system as a non-rotationally-invariant 
case, and we commented on these as we went along. To 
compute the viscosity, we shall use the strain-strain form 
of the Kubo formula, and specialize only to the ui — >• 
limit. From Eqs. (I2.9[) and p.5[) , the response function 
X is given by 

X^va^tt — > 0) = (8 va {J^PIq — Sfi,0 (Jav)o) 



lim uj 



dte iu '([V(t),^(0)]) o , 
(5.46) 



with averages taken with respect to the ground state of 
the system in the plane geometry. 

We now proceed to evaluate this expression. Introduc- 
ing momentum space creation and annihilation operators 
Ck and c£, Eq. (I5.45[) for the strain generators becomes 



J ^ - 2L 2 ^ 

kk' 



z(k-k') 



and in the limit of large system size L — > oo 
iL 2 



(5.47) 



4(2tt) 



, , d 2 kd 2 k' {k v + k' v )cl,c k x 



x (^-4) i(k " k,) ' (M8) 

(This clearly has the form of a strain generator in mo- 
mentum space.) 

To proceed further, we use the Bogoliubov transforma- 
tion 



c k = «k4 - u* k a_k, 



where 



|a k ,a k '} 



kk' 



(5.49) 



(5.50) 



and other anticommutators vanish. We have a gauge 
freedom in choosing the phases of «k and i>k, and for 
convenience we shall work in the gauge where itk is real. 
This implies that i>k transforms as an ^-wave under rota- 
tions. The explicit forms of itk and «k will not be needed. 
The Hamiltonian then takes form (see e.g. Ref. I20I ) 



Ho = E 



^ £kQ! k ak, 

k 



(5.51) 



2G 



where Eq is the ground state energy, and 



2m 



(5.52) 



is the quasiparticle dispersion relation. We work at pa- 
rameters for which £k > (gapped) at all k. 

We turn first to the contact term in Eq. (I3.5|) . Using 
the properties of the «k operators, we find 



Jnv)o — 



iL 2 



4(2tt) 



dvl 



d 2 kk v w k -i - vt-j- 



dk u 



: dk„ 



(5.53) 



This is the continuum limit of the expression given in 
Ref. H From rotational covariance, it follows that this 
term is purely antisymmetric. Writing 



«k = \Vk e 



(5.54) 



we have 



L 2 



dk. 



4(2tt 



vanish (nor does the usual inverse compressibility, which 
is defined using a derivative of pressure with volume at 
fixed particle number, not fixed chemical potential). We 
have some difficulty with the formal calculation, because 
k space was most convenient to diagonalize the Hamil- 
tonian, but to make sense of the size dependence of the 
energy, we need the formalism of Appendix with a 
confining potential (the use of periodic boundary condi- 
tions does not fit with the strain generators, though we 
could use the approach of Appendix IB]) ). But we have 
seen in previous sections how the pressure term emerges 
in the stress-stress form, and thanks to the formalism, 
this is equivalent to the strain-strain form we wanted to 
use. Hence we will not pursue this further. 

Finally, let us apply the viscosity-conductivity rela- 
tion (|4.6I) . which is still valid in the present situation, 
even though number is not conserved. We consider only 
the I = 1, or p x + ip-y, spinless superconductor. The con- 
ductivity in that case was calculated by Lutchyn et a/j*l. 
The result for the Hall conductivity (the antisymmetric 
part) in the 2D limit to order q 2 reads [cf. Eqs. (103) and 
(72) in that paper]: 



4(2^)^/ d2fck|2 ( fc ^ 

L ^ [ kdk \v k \ 2 (0 k (0 = 2vr) - MO = 0)) 



But since i>k transforms as an £-pole under rotations, we 
have 

<£ k (0 = 2vr) - k (0 = 0) = 2nl, 



a H (q,io)=I(uj) 



2 9 9 

e v F q 
4^2^+2' 



(5.58) 



where vp is the Fermi velocity and I(oS) is a dimensionless 
factor obeying J(Q)=1. Substituting this result into the 
viscosity-conductivity relation (|4.1ip . taking uj to zero, 
and remembering that the particle number density is re- 
lated to the Fermi wavevector kp by n = k F /(2n) in 
2D, we find that the Hall viscosity is i] H = hns/2, with 
s = 1/2, in accordance with the discussion above, except 
for an apparent sign discrepancy: for p + ip, s should be 
-1/2, not +1/2. 



and hence 



Jhv)q — 



L 2 



4(270"" 
Wo* 



-L 2 nse^ v . 



kdk \vk\ 



(5.55) 

(5.56) 
(5.57) 



Thus, the contact term gives the expected Hall viscosity 
r] H = hsn. 

Next, we must consider the time integral term in the 
response function X. Because of the non-degenerate 
ground state, and the gap in the excitation spectrum, 
we expect that at low frequency, the only terms will be 
the pressure and k _1 terms that go as 1/lu + . But 
here will be the derivative of the pressure with respect to 
size at fixed chemical potential fi, not fixed number N, 
and will vanish. The pressure itself is minus the deriva- 
tive of the ground state energy with volume, at fixed 
chemical potential, as usual (this "energy" is really the 
grand thermodynamic potential), and is not expected to 



VI. HALL VISCOSITY FROM 
ELECTRODYNAMICS OF A FLUID WITH 
ORBITAL SPIN 

In this section we present an alternative derivation 
of Eqs. (|1.7[) - (jl.8[> by examining the electrodynamics at 
small q and uj of a Galilean- and rotationally-invariant 
fluid in the spirit of a low-energy effective description, in 
which we assume there is an orbital spin — a per particle; 
we neglect bulk and shear viscosity. We find the orbital 
spin contribution to the conductivity, and hence to the 
viscosity. 

When one applies an electric field to such a fluid, 
to lowest order in q the electric current response is 
(j Q (q,w))(°) = a ( ^(uj)Efj(q,uj); by the continuity equa- 
tion, the corresponding change in the particle number 
density is S(n(q, uj))^ — q a (j Q (q, / 'uj + . Since 
each particle carries orbital spin — s, this leads to 
a change in the magnetization density, £/x(q, uS) = 
-si(n(q,w))' 1 '/(2m). Finally, this feeds into the elec- 

(2) 

trie current, (j„(q, w))m = tv^iq^S^ci, w), and gives the 
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following order q 2 contribution to the conductivity: 



T (2),M f 

7 vj3 



■ s q a Qn 
-% — — e 



JO) 



2 mcu 



+ C "C U a/3 



(«). (6.1) 



This contribution was previously discussed by Lutchyn 
et al<21 in the zero magnetic field case. 

A similar contribution comes from the fact that 
by Maxwell's equations, a curl of the electric field 
implies a time-dependent magnetic field, B(q,uj) = 
tupQaEfj (q, uj)/uj. The magnetic field couples to the mag- 
netization density, which is —s times the particle num- 
ber density, divided by 2m; hence, it is equivalent to an 
electric potential V B (q,uj) — ~s~B(q,w)/(2m), giving rise 
to an effective electric field E B (q,uj) = —iq tl V B (q,uj), 

(2) 

and thus to an electric current density (j„(q, w)) B = 
a^}(ui)E B (q,io). The corresponding order q 2 contribu- 
tion to the conductivity is: 

Summing Eqs. (|6.1|) - (|6.2[) we find the total contribu- 
tion of the orbital spin to the conductivity: 



J2),s 



(q,w) 



2 muj + 



JO) 



7 (0) 



where a'fhuj) is given by Eq. (|4.9|) . 



(6.3) 



These terms, however, do not in general give the full 
conductivity tensor, even for gapped systems. In that 
case, the only other possible contribution comes from a 
response with the form of the inverse (internal) compress- 
ibility k" 1 , which is not necessarily the total one as 
we will see. As before, in the presence of a nonuniform 
electric field, to lowest order in q the electric current 
response is (j a (<l,u)}^ = a^(u)Ep(q,uj). Again, the 
corresponding change in the particle number density is 
<5(n(q, oj)}^ = <Za(ja(q, uj))^/u} + . Through the inverse 
compressibility, this leads to a contribution to the pres- 
sure <5-P(q, lo) — ii~ 1 S(n(q,u})}^ /n, which is equivalent 
to an electric field E^(q,uj) = —iq fJl 8P(q,uj)/n. Finally, 
this leads to a change in electric current, (j^)jy(q, w) = 
(oj)E^(q, ui). Combining all these expressions, we 
find an order q 2 contribution to the conductivity, 



T (2),« 
T vj3 



_M"(o), 



^HflH- (6-4) 



We then have the total q 2 conductivity 

ag(q,c) = a^(q,u>)+a^ K ( q ,u). (6.5) 

Using Eq. (14.91) we then obtain, in the absence of a mag- 
netic field, 



(2)/ \ 1 



q v qp ns q" 



i 2 uj +3 2m 2 co +2 



(6.6) 



The contributions of s and k _1 have different tensor 
structures, and the result agrees with eq. (|4.15|) if we 
identify k -1 = and n H = ^ns. 

On the other hand, for nonzero magnetic field we find, 



(2),s, s ms q Q q M ns q~ 

<* v s (q,w) ^— — tape?* - tt-o— 2**13, (6.7) 

p m z u) + w c 2m z u z 

(2),/t/ s 1 qotQti ~-l 1 2 

(6.8) 

as oj — >• 0, up to terms linear or higher order in to. Adding 
these, we see that the total inverse internal compressibil- 
ity is KjlJ = +Wsuj c , and hence 



7 (2) 



(q^) 



q a q^ 



K int) 



2 ^ 1 e "P- 



(6.9) 



This agrees with the general form eq. (|4.16p , giving again 

IT 1 

T) = ^nS- 

In this derivation of the Hall viscosity, the centerpiece 
was the use of the magnetization equal to the orbital spin 
over twice the mass, for each particle. This is the stan- 
dard result for zero magnetic field, and can be obtained 
by taking minus the derivative of the kinetic energy with 
respect to B in the symmetric gauge, then setting B to 
zero. But doing the same in a uniform non-zero field 
gives —(b U 1 + h)/m for the ith particle, which is twice 
the orbital angular momentum of the cyclotron motion 
of a particle divided by the mass in the non-interacting 
case. Of course, in the rigorous microscopic derivation, 
this was not the way that the conductivity was calcu- 
lated, so there is no contradiction, however, the under- 
lying reason for the assumed coupling seems less clear. 
Hence in the magnetic field case, the present argument 
should perhaps be viewed as an interpretation of the cor- 
rect result, rather than as an alternative derivation from 
first principles. 



VII. CONCLUSION 

The motivation for and conclusions of this work have 
been described in detail in the Introduction. Essentially, 
we have a fairly complete formalism for defining the stress 
response to an external strain field, and obtaining from 
it the inverse compressibility and the viscosity tensor, 
in different situations that include non-zero temperature 
or magnetic field, and finite systems with a confining 
potential or periodic boundary conditions. The results 
from distinct formulations (such as different boundary 
conditions) agree in the thermodynamic limit. The re- 
sults illuminate the relation of the Hall viscosity results, 
obtained from adiabatic transport, to more traditional 
Kubo formula approaches; the Kubo formulas extend to 
more general situations. 
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A motivation for this work was to make it possible to 
study Hall viscosity experimentally. The generality of the 
Kubo formulas should make it possible to analyze possi- 
ble experimental set-ups to see when Hall viscosity mani- 
fests itself. In particular, the use of optical techniques to 
obtain the current response cr(q, oj) to an external per- 
turbing electromagnetic field at small q and w, in con- 
junction with the relation due to Hoyos and Son^, which 
holds more broadly as shown in Sec. HV] to obtain the 
Hall viscosity should be explored. The result would be 
of interest in quantum Hall physics, because Hall viscos- 
ity is related to the so-called shift in the ground stated, 
which can distinguish between distinct possible ground 
states at a given filling factor. 



of the magnetic field in Section (111.21) to quantize the 
Z^ v . We introduce operators Z^ and their conjugate 



momenta M.^ such that 



and take for the strain generators 

These satisfy [J Ml/ , (Zx l ) Q ] = 0, and so also 



■V,*7 (x) =0. 



(A3) 



(A4) 



(A5) 
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Appendix A: Strain generators for systems in a 
confining potential 

Because the volume of a system in the infinite plane 
with zero magnetic field is poorly defined, we would like 
to work in a finite-sized system when using the exten- 
sive forms of the stress response function given in Section 
IIII Unfortunately, because the strain generators contain 
an explicit x operator, we cannot do this with periodic 
boundary conditions. Instead, let us add to the Hamil- 
tonian H\(t) a general confining potential (this differs 
slightly from what was mentioned in Sec. III. 11 for rea- 
sons that should become clear) 



(Al) 



with an invertible matrix of shape parameters Z^, and 
we take the single-particle potential u such that u(x) — » 
oo as x — y oo in any direction. For example, the most 
general harmonic confining potential can be written 



mC 



(A2) 



Now, under an additional strain transformation we 
have x 4 — > A' T x 4 , and this A' can be absorbed into A 
by multiplication on the left. According to the definition 
just given, [/( x ) in the x variables should be unchanged, 
so we need Z — > ZA' T 1 . To implement this transforma- 
tion using operators, we are motivated by our treatment 



For given eigenvalues Z oi Z, the confining potential is 
fixed, independent of A, and defines a box that is always 
fixed in the x variables. 

To use this formalism, we first note that M. does not 
appear in the Hamiltonian H\ + U^ x \ and so with A = / 
(for example) , states with given eigenvalues Z of Z evolve 
with Z fixed, however (as with B), formal eigenstates 
of Z are not normalizable. One can, however, consider 
normalizable states that are eigenstates of Hq + [/( x ) for 
each Z, with a narrow range of Z. These are almost as 
good as true eigenstates for most purposes. 

We want to show that both our basic relation (12.21)) 
and our Kubo formulas p.3M3.5[) continue to hold unmod- 
ified in the presence of the confining potential. First, we 
will verify that the stress tensor is given by minus the 
time derivative of J^. This is simple, as J^ v commutes 
with CA X ), and the {A4,Z} terms in J Mi , commute with 



Hq X K Hence 



T$ = -i[H^+U^,J, 



(A6) 



Next, we will verify that our Kubo formulas still hold 
with these strain generators. We have the Hamiltonian 
H\(t) + U^ in the x variables (which include Z and A4), 
and we make a time-dependent canonical transformation 
to X variables. The analogs of Z and M. in X variables 
will be denoted by Z = SZS- 1 = ZA T - 1 and M = 
SAiS^ 1 . Then the Hamiltonian in X variables is H = 

jy(X) + [/(X) +Hu where 

< ) + t/ (x) = E^ + ^E^( XJ - XJ ) 



+ J2 u (-^ xi 

i 

rr _ 9 V r 
tl\ — 



(A7) 
(A8) 



to order A. In these variables, the {AA,Z} terms in Hi 
cause Z to evolve in time, in such a way that if the X*s 
also evolve by the Hi term only, then they continue to 
lie inside the "box" (defined by U) if they do so initially; 
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this was the desired behavior. The stress in these vari- 
ables, T^v , is the same expression as in the case with no 
potential U. Thus, the stress-strain Kubo formula (|3.3p 
obtained from standard linear response theory takes the 
same form as before (though the strain generator is now 
different, and the Hamiltonian includes U). Addition- 
ally, the Ward identity Eq. (|A6I) ensures that the stress- 
stress and strain-strain Kubo formulas are still valid in 
the presence of a confining potential, as is the final for- 
mula (|3.15p for the viscosity, up to a choice of the vol- 
ume L d to assign to the system. (Technical justification 
of these statements is discussed further in Appendix [C]) 
If the potential u is taken to have a hard-wall form, the 
volume can be taken as that within the walls. 

As an example of this formalism, let us consider a sys- 
tem of noninteracting spinless fermions in a harmonic 
potential in two dimensions, with Hamiltonian 



harmonically-trapped non-interacting fermions is 



H 



E 



2m 



mC< 



(A9) 



We consider the state |0) in which the potential has an- 
gular frequency C/L, meaning that |0) = L~ 1 5 l _ ll , |0) 
(strictly, this means a narrow normalizable wavepacket 
centered at this value, as we explained above), and the 
lowest Q levels filled. This state has N = Q(Q + l)/2 
particles, and its energy is 



C 



Q 

71=1 



(A10) 



We 



The integrated stress tensor is T^ v — J2iPliPv/ r 
will use the stress-stress form of the Kubo formula, eq. 
(|3.4p . Let us examine the contact term first. This gives 

(PV(0), J Q/3 (0)]) = ^ E {(pH)o S ™ + (pH)o 

i 

= (<W<^/3 + ^vot^tip) , (An) 

where the last line follows from an application of the 
virial theorem. To evaluate the time-integral term, note 
that the stress tensor and the Hamiltonian are indepen- 
dent of M^, and hence T^ v {t) leaves the ground state 
in the fixed = L~ 1 5^ v subspace. Hence, we are free 
to evaluate the time-dependence in this subspace. Doing 
so, we find 



dte^ +t {[T^{t),T afj (0)])v = 



-iE Q 1_ 

4 , _ ( u+L 
1 1 2C 



(S^aS^p + S ua 5^) . (A12) 



2lo+ 



1 - 



2C 



(A13) 

The fluid in the harmonic potential is not homoge- 
neous; its density is not uniform, and accordingly its 
other properties, such as the expectation of ^^(x), are 
not uniform either, and this effect does not disappear in 
the thermodynamic limit L, N —> oo, if we take it with 
the density at the center held fixed. A macroscopic frac- 
tion of particles experience a non-zero potential (this is 
already apparent from the use of the virial theorem: the 
kinetic energy, which is used to obtain the trace of T^, 
is only half the total energy) . Therefore use of the results 
involving P and k" 1 which relied on homogeneity is not 
justified. However, we find that if we use the free Fermi 
gas results in Sec. IV.ll to obtain the pressure and viscosity 
at a given density, and then average the results using the 
density profile of the harmonically-trapped gas, the re- 
sults agree in all details with the thermodynamic limit of 
the above. To obtain a fluid that is homogeneous (up to 
a negligible boundary layer) in the thermodynamic limit, 
we need a confining potential that is essentially zero in 
the interior, then rises rapidly very close to the edge. 



Appendix B: Stress response without use of strain 
generators 

In this section we show how to obtain even more gen- 
eral forms of the stress response, which work for either 
the infinite system with a confining potential, or for pe- 
riodic boundary conditions, and are analogous to the 
stress-strain and strain-strain forms, but do not require 
the use or existence of the strain generators. The ap- 
proach we use is based on that of Niu et al. for the 
conductivity case22, however we extend it to non-zero 
frequency as well as adapting it for the stress response. 
For the case of infinite space, the expressions reduce to 
the forms in the text when written in terms of strain 
generators. 

We begin with the time-integral part of the extensive 
stress-stress form at zero temperature, which we write 
in a spectral representation. For definiteness, one can 
assume periodic boundary conditions. Then we have 

1 r°° 



Hence we find that the response function for 



dte^ t {[T lw {t),T aP ^)]) Q = 

1 r°° r 
— / dte^ +t V \ e -^- E ^ (0|T^(0)|n)(n|r Q ^(0)|0) 
^ Jo ~ 

- e-^°- B «>* (0\T aP {0)\n)(n\T^(0)\0)] , ( B1 ) 

where {\n)} is an orthonormal set of energy eigenstates 
of H (or H + U if the confining potential is used in place 
of periodic boundary conditions) with energies E n , and 



30 



we now write |0) for the ground state, which we assume 
for simplicity is non-degenerate. We now use the identity 



-i(E n -E )t 



d_ -i(E„-E )t 
d£_ 

-i(E n - E ) 



(B2) 



in the first term inside the integral, and the same with 
n and switched in the second term; then integrate by 
parts. Further, we recall that in the x variables, the 
Hamiltonian H\(t) to first order in A is 



H A {t)=H o -T a0 X a0 {t) 



(B3) 



(where T a p = T a( g(0) is at zero time in the Heisenberg 
picture), and so if we define the ground state of H\ for 
given A to be \ip(X)), that is 



H A \p(X)) = E(X)\cp(X)} 



(B4) 



(with \<p(0)) = |0), E(0) = E ), then perturbation theory 
to first order in A gives us 



dp 



dX af: 



E 



|n)(n|r^(0)|0) 
E n — E 



(B5) 



(We leave it as understood that the partial derivative 
is evaluated at A = 0.) Combining these results, and 
moving the term resulting from the lower limit in the 
integration by parts to the other side, we arrive at two 
forms of the response function X^ a p{uj): 



<P(P) 



V(0) 



dtp 
dX~ 



a/3 



dip 



a0 



<p(0) 



— dte 
o 



dt e 



- ht M0)|[T M!/ (t),T^(0)]|^(0)) (B6) 
dip 



<P(0) 



dp 



dX a /3 



<P(P) 



(B7) 



,dX a p 

The first of these is the stress-stress form, and the second 
is the stress-strain form. 

In order to use the same identity and perform a second 
integration by parts, we generalize the previous pertur- 
bation formula by using T uv (t) for given t as the pertur- 
bation, and define 

Mt)\ = y \n)(n\T^(t)\0) 



E 



E n — E 
\n)(n\T^(0)\0) i(Eo _ En)t 



E n - E 

Then the strain-strain form of X is 

X uva f}(uj) = 



(B8) 
(B9) 



The first part is the curvature of the Berry connection. 
This part reproduces the formulas for the antisymmet- 
ric (Hall viscosity) part of X at ui = that result from 
adiabatic transport in a gapped system (with the sign 
here corrected as explained earlier). The derivation of 
this w = limit from the stress-strain form is presum- 
ably equivalent to the one reviewed in Ref. |5|. The last 
term produces the in /cj^ term in the trace part, as 
uj —> 0, similarly as in Sec. IV. 21 (such parts, referred to 
as "persistent currents" by analogy with the conductiv- 
ity case, were suppressed in Refs. [a by subtracting the 
ground state energy E(X) from H\). The strain-strain 
form given here can be viewed as generalizing the adia- 
batic curvature (and the related Chern number for con- 
ductivity) to the full stress response tensor, to all fre- 
quencies, and to systems that are gapless in the thermo- 
dynamic limit. Because the formal derivation is general 
in form, similar ones can be given for conductivity at all 
u>, and for other transport properties also. 

The similarity with the three forms of X given in the 
text is evident, and can be made exact by use of the 
formula 



- OJ 



dt e 



(BIO) 

MPi\_/MPi d< p(°) V 
dX aP I \ dX a p dX aw I 

( /dtp(t) cV(0)\ /dp(0) dip(t)\ 



( d<p(0) 



\ dX^ dX a p / \ dX a/3 dX^ J 



lv(o)), 



(Bll) 



and expanding to first order. Then one can view 



= -iJ aP (t)\v(0)) = 

dX a p I 



E 



\n)(n\T afi (t)\0) 
E n — Eq 



(B12) 

as resulting from the basic relation, eq. (|2.21j) . This is 
not merely an analogy; for the infinite-space geometry, 
the ground states of H\ are of exactly this form (with 
t = 0), using the formalism for magnetic field and for a 
confining potential, so they form a "homogeneous bun- 
dle" in the language of Ref. [H (here for GL(<i, R), not just 
SL(d, K) as mainly considered there). We note however 
that strictly speaking for the infinite-space geometry, the 
formalism for a magnetic field or in Appendix [X] again 
produces issues regarding the existence of normalizable 
energy eigenstates, and some arguments with wavepack- 
ets analogous to those in the text or in Appendix [C] are 
needed to overcome these. 

In any case, the forms given here for periodic boundary 
conditions are expected to yield the same results for in- 
tensive quantities in the thermodynamic limit, justifying 
the use of the same notation X nva p(uj) in the formulas 
given here. 



Appendix C: Time-translation invariance and the 
Kubo formulas 



In this Appendix, we give arguments that justify the 
use of time-translation invariance (TTI) in the derivation 
of the Kubo formulas in Sec. IIII.ll even though the states 
that must be used are not true eigenstates of the Hamil- 
tonian. We say that a correlation or response function in 
the time domain has TTI if shifting the time argument 
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in each operator by the same constant has no effect on it. 
For the zero magnetic field case in the informal treatment 
without a confining potential that we have mentioned in 
Sec. III. 1 1 one has to appeal to an assumed limit in which 
time-dependence of the state occurs only near the bound- 
ary, and the bulk of the system presumably dominates 
the response over relevant time scales. This is intuitively 
appealing, but quite involved to justify fully. Here we 
give more detailed arguments for the case with a mag- 
netic field, in which normalizable ground (energy eigen-) 
states in which a disk-shaped region is occupied exist for 
given magnetic field. There is still an issue here, however, 
because our strain generators involved the introduction 
of the magnetic field variable B, and eigenstates of B are 
not normalizable. Similar arguments are also given for 
the case with a confining potential, in the formalism of 
Appendix [2] 

The "ground" state that we described in Sec. 111.21 is 
not strictly an energy eigenstate, because though it is 
such for each value of B, the energy eigenvalue depends 
on B (the Hamiltonian Hq depends on B but not on V, 
so parts of the state with different B values do not mix) . 
Thus the contributions from different B change phase 
with time at different rates. When we calculate an ex- 
pectation value of some operators, these time-dependent 
phase may cancel; in particular, they will if those oper- 
ators do not contain V . Now the stress tensor T^ v does 
not contain V, so the correlation function 

([T^(t),T a p(t')}) =f(t,t') (CI) 

is TTI — it is unchanged if we replace t, t' by t + to, t' + <o- 
(In the following, the indices on the operators will play 
no role, so they will not be recorded on / and similar 
functions below.) Now consider the correlation function 

([T li u(t),J a0 (t')]) o = g(t,t'). (C2) 

From the identity (1 2 . 5 6 [) . 

JLg(t,t') = -f(t,t% (C3) 

and f(t,f) = f(t - t') (say) by TTI. Integrating gives 

g(t,t')= f 1 dt"f(t")+g(t,t), (CM) 
Jo 

where 

g(t,t) = (e iHot [1V(0), J a p(0)]e- iHot ) Q . (C5) 

The commutator [T M „(0), J Q/ 3(0)] does not contain V for 
any choice of indices, and so g{t,t) is again independent 



of t, or TTI. Hence g(t, t') = g(t - t'), say, is TTI. These 
results justify the use of TTI to pass from the stress- 
strain to the stress-stress form of response. 

For the strain-strain form, we proceed similarly with 

([•Mi), J a p(t)]) = h(t,t'). (C6) 

Then 

^-h(t,t') = -g(t-t'), (C7) 

and so 

pt-t' 

h(t, t')= - I dt" g{t") + h(t, t). (C8) 
Jo 

The equal-time piece is 

h(t,t) = ([J liU (t),J a0 (t)]) o . (C9) 

The commutators of two J's are the gl(d, R) Lie algebra 
relations, given in (|2.9[) . The generator corresponding to 
the trace (the generator of the gl(l) or u(l) subalgebra) 
never occurs on the right-hand side; see eq. (|3.29p . The 
traceless parts of J's do not contain V, and so h(t,t) 
is independent of t, that is, it is TTI. Hence h(t,t') = 
h(t-t'), say is TTI. 

In fact, TTI was not used to obtain the strain-strain 
form from the stress-strain form, but we have learned 
that it can be applied to the final form. This enables us 
to reverse the argument, but with the roles of the two 
pairs of indices interchanged. This leads to the same 
stress-stress form of response, but with the two pairs of 
indices exchanged in the contact term. We showed in Sec. 
IIII.5I that this contact term is symmetric, which shows 
the arguments are consistent. 

For the formalism of Appendix [5] with a confining po- 
tential, we can proceed similarly, using a wavepacket in 
Z space of energy eigenstates for each Z. The preced- 
ing line of argument goes through, up to the point where 
h(t, t) was studied. In the present case, all strain gener- 
ators contain AA, not only the trace. Hence h(t,t) is 
not independent of t. If we take its time derivative, we 
find 

^h(t,t) = iS au (T^(t)) - i6^(T al/ (t)) , (CIO) 

in which the trace of the expectation of the stress can- 
cels. We have shown that the expectation of the traceless 
part of T^u is time independent, so h(t,t) is linear in t, 
and also its time derivative is smaller than 0(L d ) as the 
thermodynamic limit is taken. Thus in the limit, for our 
purposes h(t,t') is again TTI. The strain-strain form in 
the main text is nonetheless correct as written even in 
finite size. 
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